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ABSTRACT 

We have obtained spectra of 135 H II regions located in the inner and extended disks of the spiral galaxies 
NGC 1512 and NGC 3621, spanning the range of galactocentric distances 0.2 - 2 x R25 (from ^2-3 kpc to ^18- 
25 kpc). We find that the excitation properties of nebulae in the outer {R > R2;,) disks are similar to those of the 
inner disks, but on average younger H II regions tend to be selected in the bright inner disks. Reddening by dust 
is not negligible in the outer disks, and subject to significant large-scale spatial variations. For both galaxies 
the radial abundance gradient flattens to a constant value outside of the isophotal radius. The outer disk O/H 
abundance ratio is highly homogeneous, with a scatter of only ^0.06 dex. In the case of the interacting galaxy 
NGC 1512 we find a number of Hll regions with peculiar metallicity for their radius, a result which can be 
interpreted by gas flows activated by the gravitational encounter with NGC 1510. Based on the excitation and 
chemical (N/O ratio) analysis we find no compelling evidence for variations in the upper initial mass function 
of ionizing clusters of extended disks. The O/H abundance in the outer disks of the target galaxies corresponds 
to ~35% of the solar value (or higher, depending on the metallicity diagnostic). This agrees with our earlier 
measurements in M83 and NGC 4625, and conflicts with the notion that metallicities in extended disks of spiral 
galaxies are low and on the order of ~'0. IxZq. We show that, in general, the observed metal enrichment cannot 
be produced with the current level of star formation, even if the latter extends over a Hubble time. We discuss 
the possibility that metal transport mechanisms from the inner disks lead to metal pollution of the outer disks. 
Gas accretion from the intergalactic medium, enriched by outflows, offers an alternative solution, justified 
within the framework of hydrodynamic simulations of galaxy evolution. Specific model predictions of the 
chemical enrichment and the flat gradients in extended disks of nearby galaxies will be valuable to discriminate 
between these different scenarios. 
Subject headings: galaxies; abundances — galaxies; ISM — galaxies; individual (NGC 1512, NGC 3621) 



1. INTRODUCTION 

Optical and far-UV observations reveal that a large fraction 
of star forming galaxies in the local Universe harbor signifi- 
cant levels of recent star formation in their exten ded disks, in 
the form of young star clusters (ages < 1 Gyr; Dong et al. 



2008 



Herbert-Fort et al.||2009|l distributed well beyond the 



optical edges of the disks (Gil de Paz et al."2005'i 'Zaritsky' 
& Christlein 2007; Thilker et al. 2005, 2007; Lemonias et al. 
201 1 ). This disco very h as prompted several recent studies 
Te.g. .Boissier et aL][2n07] iBigiel et al..2010b^,Goddard et al., 



2010 ' Alberts et al.||20lib concerning the conversion of gas 
mto stars in a poorly explored galactic environment, charac- 
terized by low gas densities, long dynamical timescales, and 
presumably a lower chemical enrichment compared to the in- 
ner disks of galaxies. Remarkably, also a number of early- 
type galaxies, which are generally quiescent in terms of star 
formation, have b een found to host UV-emitting star clusters 
in flieir outskirts (^aJim & Rich||2010l [Thilker et al.||2010t 
|Moffett et al.|201 l| l, and thus appear to be experiencing a 're- 

' Based on observations collected at the European Southern Observa- 
tory, Chile, under program 084.B-0107. 



juvenation' in their extended disks. 

A variety of mechanisms can be invoked for the assembly 
of the large H I envelopes observed around outer disk host 
galaxies, providing the raw material for the formation of new 
stars, and restructure the gas distribution into the observed ex- 
tended, star forming disks. Among these are tidal interactions, 
mergers, bar-induced gas flows and gas accretion from the in- 
tergalactic medium. In particula r, 3/4 of the Type 1 UV disks 
defined by Thilker et al. ( 2007) l display signatures of recent 
or ongoing gravitational encounters. However, the processes 
that trigger and sustain the star forming activity in outer disks 
are not always obvious, and still poorly understood in the case 
of isolated galaxies. Once the gas reservoir is in place, spi- 
ral density waves are though t to be able to propagate from 
the inner into the outer disks ( |Bertin & Amorisco|2010| ) and 
trigger star formation in overdensities of the gas distribution 
(Bush et al. 200 8, 2010), even if the spatially-averaged gas 
density r emains b elow the critical value for gravitational in- 
stability ( Dong etal. ,2008) . However, this is unlikely to be the 
only process involved in the generation of extended UV disks, 
even for isolated galaxies, becau se the phenomenon is also 
observed around dwarf galaxies ( [Hunter et al.||20lT] l, where 
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spiral arms are generally absent. 

The star forming activity in extended disks proceeds with 
a very small efficiency, with values of the star formation 
rate surface density that are one to two orders of magnitude 
smaller than those observed in the inner disks, with deple- 
tion timescal es that are on the order of 10" years (Bigielj 



et al. 2010b| l. Under these conditions a very slow buildup 
of metals can be reasonably anticipated, but the present-day 
oxygen abundances measured in outer disks seem to contrast 
this expectation. In the few detailed studies of the chemi- 
cal composition of the interstellar medium (ISM) of extended 
disks carried out so far from H II region spectroscopy, typical 
O/H values corre sponding to at least 0.3-0.4 times the solar 
valu^are found ( |Bresolin et al.|2009bl [Goddard et al.|20I I 
Werketal.|20IIi), which are difficult to reconcile with the sig 



niticantly lower metallicities predicted by simple models of 
galactic evolution for systems that, like extended disks, have 
extremely high gas fractions (> 90%). In addition, these pio- 
neering studies have also found that beyond the isophotal radii 
(7?25) the metallicity radial distribution flattens to a virtually 
constant value, irrespective of the galactocentric distance, in 
contrast with the negative exponential radial abundance gradi- 
ents present at smaller distances. This behavior has also been 
seen in the old stellar populations ( Vlajic et al. 2009). 

Various processes of angular momentum redistribution and 
gas mixing in galaxy disks have been invoked to explain these 
observations, including resonance scattering with transient 
spirals (Roskar et al. 2008b) and spiral-bar resonance over- 
lap (Minchev & Famaey 2010). Alternatively, Salim & Rich 
pOlO) and Lemonias et al.^(2011) , among others, have argued 
that the gas that provides the raw material for the ongoing star 
formation in extended disks results from cold-mode accretion, 
as predicted by cosmological simulation s (Keres et al. 2005; 
Dekel & Bimboim|2006||Ocvirk et al^DOSp . This idea, how- 
ever, contrasts with the low metal content expected for the 
inf ailing gas. In order to explain the relatively high metal- 
licity found for the ISM of the outer disks we would also 
need to consider a mechanism similar to the enriched infall 
scenario predicted b y simulat ions involving win d recycling 
( [Oppenheimer & Da ve 2008; Dave et al.|2"0TTb] l. We can of 
course envision a combination of internally- and externally- 
driven processes to enrich and homogenize the ISM at large 
galactocentric distances. 

Measuring chemical abundances is of primary importance 
in the study of outer disks, in order to provide observational 
tests for the mechanisms mentioned above, and to understand 
whether extended disks truly represent young, chemically un- 
evolved portions of galaxies, as implied by the inside-out disk 



growth para digm (Matteucci & Francois 1989[ Boissier & 
|Prantzosp9 99 ; Naab & Ostriker 2006 ), and whether secular 
processes, too, affect the radial distribution of metals in the 
external regions of galaxies. 

While only the youngest stellar clusters, still contain- 
ing short-lived massive ionizing stars, produce a measurable 
emission line spectrum, the optical spectroscopy of the associ- 
ated ionized gas provides a unique probe into the present-day 
metallicity of outer galactic disks. In this paper we continue 
our exploration of the chemical abundances of the faint H II 
regions located in the extended disks of nearby spiral galax- 
ies. In Bresolin et al. (2009b) and Goddard et al. (2011 ) we 
presented a chemical abundance analysis of the prototypical 



^ We adopt the solar oxygen abundance 12 + log(O/H)0 =8.69 bylAs- 
undetal.|p009) . 



outer disk galaxies M83 and NGC 4625, respectively. Here 
we apply the same multi-object spectroscopic technique to 
obtain the chemical abundances of a large number of H II re- 
gions in the outer disks of the southern hemisphere galaxies 
NGC 1512 and NGC 3621. Our radial coverage extends from 
the inner disk (--0.2 x R25) out to ^2 x R25 (19-25 kpc) for 
both galaxies. The comparison of their chemical properties is 
of particular interest, because while NGC 3621 is a bulgeless, 
isolated spiral, NGC 1512 represents, together with its com- 
panion NGC 1510, an interacting system. Therefore, while 
both galaxies display a multitude of UV-bright clusters and 
H II regions in their extended disks, their origin and chemical 
properties can in principle differ, since different mixing pro- 
cesses could be at work. Further details about the two galax- 
ies are introduced in Section 2, together with a description of 
our observations and the data reduction. In Section 3 we pro- 
vide an account of the excitation properties of our H II region 
sample. The chemical abundances and the radial abundance 
gradients are derived in Section 4. Our interpretation of the re- 
sults obtained in this work, and more generally in our project 
on extended disks, in terms of metal mixing and enriched ac- 
cretion, follows in Section 5. We summarize our results and 
conclusions in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

In this section we first provide a brief description of the 
properties of the two galaxies investigated in this paper that 
are relevant for our analysis. The observations and the data 
reduction are presented next. 

2.1. The target galaxies 

NGC 1512 - This galaxy and its companion, the blue 
compact dwarf NGC 1510 (5' or -13.8 kpc to the SW), are 
a well-known interacting pair The tidal interaction is likely 
responsible for the nuclear hotspot ring, the bar structure, the 
larger inner ring knotted with H II regions (hence the SB(r) 
morphological classification), and the star forming activity 
seen in the asymmetric, filamentary outer spiral arm structure, 
composed of a large number of UV-bright clu sters and Hn_ 
regions (Kinman 1978; Hawar den et al.|[T979[ |Thilker et aE] 
2007; see Fig.fTll. The ongoing star formation m the compact 
NGC 1510 has also likely been triggered by the interaction. 
Evidence for tidally stripped material can be found in the 
zone between the two galaxies and in the wispy filaments 
to the NW of the pair Koriba lski & L6pez-Sanchez| p009| l 
made a detail ed investigation of the e xtended Hi envelope 
discovered by [Hawarden et al.| ( |1979| l, and discovered two 
tidal dwarf galaxy candidates at galactocentric distances 
between 64 kpc and 80 kpc. The spatial extent of the H I 
distribution (—55 kpc in radius) is comparable to what is 
seen in other extended disk galaxies, such as M83. T he hy- 
drodynamical simulations of the galaxy pair by ChakrabartiJ 
poll) indicate a mass ratio of 1:50. A dynamical mass 
Gm5l2 of 3 X 10" M© was measured by KoribalsM] 



|P 



& Lopez-Sanchez (2009 1. These authors have also found a 
good spatial correlation between the Hi emission and the 
distribution of UV clusters in the Galaxy Evolution Explorer 
(GALEX) images, with the clusters located at the peaks of 
the gas column density. Star fo rmation occurs near the local 
critical density ( :Kennicutt||1989 ). Ha images of NGC 1512 
(and also NGC 3621) show that the spatial extent of the H II 
region distribution is similar to that of the extended UV disk. 
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Fig. 1 . — Near-UV GALEX image of NGC 1512, identifying the target H II regions. The ellipse represents the projected circle with radius equal to the isophotal 
radius R25 . The companion galaxy NGC 15 10 is identified in the south-west direction. The wedge-shaped area delimited by the blue lines marks the 'tidal bridge' 
between NGC 1510 and NGC 1512. Regions num bered from 53 to 57 are located along this feature. The inset shows the inner 3f2 of the galaxy from the 
continuum-subtracted SINGG (Meurer et al.|2006) Ha image. Region numbers correspond to those given in TablelS] 



NGC 3621 - The GALEX and Ha images (the latter is shown 

in Fig.|2| show a rather regular spiral structure extending out 

to a large galactocentric distance, and conta ined within the 

34 X 7 kpc H I envelope ( [Koribalski et al. 1200 4). This galaxy 

is less massi ve than NGC 1512. with a dynamical mass of 
.„in.. IT — T-rr . ■ „i .,„ — ■' — 



2 X 10^" Mq (deBloketal. 2008). NGC 3621 appears to be 



isolated ( Kraan-Korteweg & Tammann_1979j. While the H i 
velocity maps suggest past interactions ( |Thilker et aH2007ll, 
rotation in a warped disk offers an alternative explanation (|de| 
Blok et aLpOOS ). An analysis of GALEX and Spitzer images 



of the outer disk of this galaxy is presented by Alberts et al. 
(201 1 1, who derived cluster ages between 1 Myr and 1 Gyr. 



2.2. Observations and data reduction 

Table [1] summarizes the main properties of the target galax- 
ies usedfor this study. The center coordinates, as well as 
the positions of the individual H 11 regions, were measured 
from astro metrically calibrated Ha images obtained by the 
SINGG ( Meurer et al.||2006| l and SINGS (Kennicut t eTaT] 
[2003b ) projects for NGC 1512 and NGC 3621, respectively 
(the preimaging data described below were used to measure 
the coordinates of a few Hll regions not included in the 
SINGG Ha image of NGC 1512). For NGC 3621 we adopted 



the Cepheid distance given by'Freedma n et al?] (2001 1 , while 
for NGC 1512 the distance has been derived, following|Korib- 



TABLE 1 
Galaxy parameters 



NGC 1512 



NGC 3621 



R.A. (J2000.0) 04:03:54.17 11:18:16.52 

Decl. (J2000.0) -43:20:55.4 -32:48:50.7 

Morphological type SB(r)a SA(s)d 

Distance 9.5 Mpc 6.6 Mpc 

R25 4^26(11.76 kpc)" 4.'89 (9.39 kpc)'' 

Inclination (deg) 35"^ 65'' 

Position angle (deg) 260"^ 345'' 

B^. 10.96'= 9.20"= 

M° -18.93 -19.90 

N ote. — "1 arcse c = 46 pc. * 1 a i -csec = 3 2 pc. '^ [K^balsk^& 
L6pez -Sanch"ez]j2009| . "^ [de Bloket ar]j2008t . ' |de Vaucouleurs et al. 

^i99fr — — — 



alski et al. ( 2004| l, from the recessional velocity and a value of 
Ho = 75 km s"' Mpc"' . The position angle and inclination to 
the line of sight of NGC 1512 corre spond to the H 1 kinemat- 
ics value s determined by Koribalski & Lopez-Sanchez (2009] 
see also Hawarden et al.||1979|), w hic h agree with the opti - 
cal measurements by IButa 1 1988 1 and Dicaire et al. ( 2008[ ). 
For NGC 3621 we adop t ed the average values from the Hi 
study by |de Blok et al.| P008| l. The remaining parameters 
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Fig. 2.— Conti nuum-subtracted Ha image ofNGC 3621 (SINGS, [Kemir] 
|cutt et al.|2003b) , identifying the target H II regions. The ellipse shows the 
projected circle with radius equal to the isophotal radius R25 . Region num- 
bers correspond to those given in Table|6] 



included in Table [T] have been collected f romlde Vaucouleursl 
|etal.| (|1991) and the HyperLeda database ( [Paturel et al.|2003| l. 

We obtained optical spectroscopy of H II region candidates 
with the Focal Optical Reducer and Spectrograph (FORS2) at- 
tached to the 8 m-aperture Antu unit of the European Southern 
Observatory (ESO) Very Large Telescope on Cerro Paranal. 
The observations were carried out in service mode between 
November 2009 and March 2010. The targets for the multi- 
object spectroscopy were selected from narrow-band Ha pre- 
imaging data taken with the same instrument. The instrument 
setup was the same o ne we used in our s tudy of the outer disk 
H II regions of M83 (jBresolin et al.'2009b). Masks were cut 
with 1-arcsec wide slits. ITie 600B grism with a central wave- 
length of 4650 A provided coverage of the blue spectral range 
(~3600-5200 A, at 4.5 A FWHM spectral resolution), while 
the 1200R grism and a central wavelength of 6500 A allowed 
us to cover the region around Ha, [N ll] AA6548, 6583 and the 
[S II] AA6717, 6731 lines with a spectral resolution of 2.4 A. 

The targets are distributed in four 6.'8 x 6^8 F0RS2 fields 
in NGC 1512, and in three fields in NGC 3621. H II regions in 
both the inner and outer disks were included in the slit masks, 
to ensure a wide radial coverage for the resulting nebular 
chemical abundance gradient. In the final selection of the tar- 
get H II regions we accounted for both the Ha brightness and 
the spatial distribution, observing some of the brightest H II 



TABLE 2 
Positions of SNRs and Emission Line stars 

R.A. (J2000) DEC (J2000) Comment 

Supernova remnant candidates 

0403 46.90 -43 20 32.8 inNGC1512 

04 03 46.00 -43 22 44.4 

111809.39 -324320.9 inNGC3621 

1118 12.02 -32 46 50.0 

Emission-line stars (NGC 3621) 

1118 24.15 -32 50 09.8 

11 18 19.43 -32 56 01.4 

1118 12.00 -32 42 06.5 

Note. — Units of right ascension are hours, minutes and seconds, 
and units of declination are degrees, arcminutes and arcseconds. 

regions in the outer disks, but including several fainter ones 
in order to fill the slit maks. The exposure times adopted were 
6x940 s and 6x970 s for the blue exposures in NGC 1512 and 
NGC 3621, respectively. In the red 3 x 667 s exposures were 
used. For the easternmost field of NGC 1512 no red exposure 
could be obtained during the time allocated for this program. 
Therefore, for 10 targets in NGC 1512 the Ha, [N ll] and [S ll] 
flux information is unavailable. We also secured spectra of the 
spectrophotometric standard stars LTT 1788 and LTT 3864, 
in addition to those included in the Paranal Observatory cal- 
ibration plan, in order to determine the instrumental spectral 
response. 

The raw data were reduced with the EsoRex pipeline pro- 
vided by ESO. This first step of the data reduction yielded 
wavelength-calibrated, two-dimensional spectra for each in- 
dividual target. The final extractions and the flux calibration 
were obtained with standard iRAlQ routines. The emission 
line intensities were then measured with the non-interactive 
f itprof s program, by fitting one-dimensional gaussian pro- 
files to the most important spectral features. The line intensity 
errors were estimated by accounting for both the uncertainty 
in placing the continuum level, as provided by the rms in re- 
gions adjacent to the emission lines, and the photon statistics. 

A number of targets were removed because they were either 
too faint to provide useful line emission information, or they 
were found to be background galaxies. In addition, we ex- 
cluded from our analysis four objects (two in each galaxy) that 
we classified as supernova remnants (SNRs) from the strength 
of their [Ol] and [S ll] lines relative to Ha. In NGC 3621 
we also identified three emission-line star candidates, whose 
spectra are characterized by broad Ha lines, similar to those 
id entified in the first s tellar spectroscopic study of this galaxy 
by'Bresolin et al.|p00 1). The coordinates of these objects and 
of the SNRs are reported in Table l2] 

Our final sample comprises 62 and 73 Hll regions in 
NGC 1512 and NGC 3621, respectively. Their celestial coor- 
dinates are presented in Appendix A (Tables l5] and l6|l, where 
the targets are listed in order of decreasing declination, while 
their spatial distributions are shown in Fig.[T]and Fig.|2] 

3. EXTINCTION 

The blue and red spectra have no overlap in wavelength. 
For this reason, the Ha/H/3 ratio could not be used to derive 

' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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Fig. 3.— The c(H^) distribution across NGC 1512 (left) and NGC 3621 (rigiit). 
proportional to tlie measured extinction. 



Each H II region in our sample is represented by a circle with a radius 



the value of the interstellar reddening from the B aimer decre- 
ment, and the strongest Balmer lines present in the blue spec- 
tral ran ge (H/?, H7 an d US) were used for this purpose, adopt- 
ing the Seaton| ( [19791 ) extinction curve. We accounted for the 
effects of underlying stellar absorption, by requiring that the 
H7/H/3 and HSfHf3 ratios provided the same value for the ex- 
tinction coefficient. In the few cases (^ 6% of the sample) 
where HS was too faint to be measured reliably, we assumed 
an equivalent width of 2 A for the absorption component. 

The reddening-corrected fluxes for the [O ll] A3727, 
[Oiii]A5007, [Nii]A6583 and [Sii] AA6717,6731 emission 
lines are summarized in Tables |5] and |6l The associated er- 
rors include the line intensity uncertainties estimated above, 
and a term obtained by propagating the uncertainty in c(H/3), 
the logarithmic extinction at HpTj For the normalization 
of the [Nil] and [S ll] line fluxes to H/3= 100 we assumed 
Ha = 286, which is valid for case B emission at 10"* K. 
The two tables also report the c(H/3) values, and the to- 
tal, extinction-corrected Ha fluxes, obtained from the flux- 
calibrated narrow-band images used for the astrometry (a cor- 
rection for the [N ll] contamination to the emission measured 
by the narrow-band filters was evaluated from our spectra). 

We found that in the outer disks of the two galaxies studied 
here the extinction is not negligible (we note that the fore- 
ground reddening values according to Schlegel et al. 1998] 
are quite small, E(B-V) = 0.011 and 0.08 for NGC 1512 and 
NGC 3621, respectively). Outer disks are therefore not dust- 
free, a result that is supported by the relatively high oxy- 
gen abundance we find (Section 4). We also find that the 
extinction is variable on large spatial scales. In particular, 
in the eastern outer arm of NGC 1512 the median extinc- 
tion value is c(H/3) = 0.34, while on the opposite side of the 
galaxy it is essentially zero. A similar asymmetry is present 
in NGC 3621, where the northern outer arm appears less dusty 



than the southern one. For the latter the section extending be- 
tween regions 59 and 71 has a median c(H/3) = 0.37, peaking 
at a value around 0.50 in the arm section between regions 59 
and 65 (separated by 1.8 kpc in projection). Fig.|3]shows the 
spatial distribution of the c(H/3) values across the disks of the 
two galaxies. These asymmetrical distributions reflect those 
in the neutral gas density, as c an be seen from the Hi maps 
sho wn by [K oribalski & L6pe z-Sanchez| ( [20091 NGC 1512) 
and lWalter et al., ( |2008, NGC 352TT 

4. EXCITATION PROPERTIES 

We inspected the excitation properties of our H II region 
sample by plotting the [N ll] A6583/Hq; (=N2) fine ratio 
against [O III] A5007/H/J (=03), as shown in Fig.|4|(the left 
and right panels refer to NGC 1512 and NGC 3621, respec- 
tively). For comparison, we have included the samples of in- 
ner and oute r disk H 11 regions studied in M83 by Bresolm] 
[eTaL] ( [2009 b) and in NGC 4625 by Godda rd et al.| ( |201I| 7 
and shown with triangle symbols. Additional data points for 
nearly 500 extragalactic inner disk H 11 regions in spiral and 
irregular galaxies (cross symbols) have been drawn from th e 
following paper s: McCafl et al. (1985), Garn ett et al.| ( |1997| ), 
I BresoUn et aL| [1999| |2005| |2009a|), |van Zee et al.| P^ , 
jCennic utt et al.| ( |2003a| ), and| Esteban et al. |(| 2009 p. Lastly, 
we also included the 39 H 11 regions located in the Milky Way 
and the Magellanic Clouds studied by ,Kennicutt et al.,(2000[ 
square symbols), because in many cases these nebulae only 
contain a single or very few ionizing stars, a situation re- 



E(B-V) = 0.68 c(H^) Seaton 1979 . 



sembling the one encountered in outer disks (Bresolin et al. 
2009b). For both NGC 1512 and NGC 3621 the area oTtEe 
symbols is proportional to the observed Ha flux, in order to 
provide information on the relative ionizing output of the em- 
bedded clusters. 

As an aid in the interpretation of the excitation diagrams, 
we have fitted the sequence outlined by the data points with a 
four-degree polynomial, obtaining; 
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Fig. 4. — The [Nii]/Hq: vs. [Oiii]/H/3 diagnostic diagram for tlie Hll regions observed in NGC 1512 (left) and NGC 3621 (right), divid ed into inner (blue 
circles) and outer (red circles) disk objects. The samples of ionized nebulae analyzed in the extended disk galaxies M83 jBresolin et al.|2009bj and NGC 4625 
JGoddard et al. 2011 1 are included (triangles). Comparison samples from the Milky Way and the Magellanic Clouds (squares) and additional galaxies (crosses) 
are shown (see text tor references). The two curves enclose H II regions that deviate by ±0.15 dex from the fit representing the empirical excitation sequence. 
The sizes of the symbols shown for NGC 1512 and NGC 3621 are proportional to the observed Ha flux. 
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where x = log([N ll] A6583/Ha), y = log([0 ill] A5007/H/3). In 
Fig. HI we show two curves, separated by ±0.15 dex from this 
empirical sequence, defining a region that encloses the bulk 
of the data points. Some important deviations from the main 
sequence occur both above and below this zone. The upper 
region of the diagnostic diagram hosts nebulae with unusually 
hard ionizing fields. In fact, in virtually all cases, H II regions 
located in this area of the diagram display, in addition to large 
[O lll]/H/3 flux ratios, also nebular He II A4686 emission, as 
summarized recently by Bresolin ( 2011b ). We note that for 
our target 1 in NGC 3621, identified in Fig. |4] (right) as the 
outlier located in the central top portion of the diagram, we 
also have detected the He II A4686 emission line, and thus this 
nebula represents a new instance of this rare class of high- 
excitation extragaiactic H II regions, whose ionizing sources 
are still a matter of debate ( [Kehrig et al.|201 1; ). (The fact that 
we can spatially resolve this target in our Ha images excludes 
the possibility that this high-excitation nebula is a planetary 

nebula)^ 

Goddard et al. ( 2011| l noticed that a significant fraction of 
outer disk H II regions in NGC 4625 (green triangles in Fig.|4| 
lies below the empirical excitation sequence in the N2 vs. 03 
diagram defined above. Fig.Hlshows that this region of the di- 
agram, characterized by a lower excitation ([Olll]/H/3 ratio) 
compared to the main sequence, is populated also by inner 
disk H II regions (cross symbols), even though these cases are 
apparently rare. The few Hll regions in NGC 1512 that He 



below the main excitation sequence belong to both the inner 
disk (defined as R/R25 < 1, and represented by blue symbols) 
and the outer disk {R/R25 > 1, red symbols) of the galaxy, and 
tend to have smaller ionizing fluxes compared to the major- 
ity of the objects located along the sequence, as indicated by 
the size of the symbols. This situation resembles the case of 
the H II regions in the outer disk of M83 (orange triangles), 
where only very few points deviate from the main excitation 
sequence. 

The case of NGC 3621 is quite different. About half of 
the Hll regions in this galaxy, in both the inner and outer 
parts of t he disk, deviate sign ificantly from the excitation 
sequence. Goddard et al. (2011) attributed the low excitation 
of outer dlisk Hll regions in NGC 4625 to a combination 
of two effects. First, the progressive depletion of massive 
stars belonging to the upper main sequence of the ionizing 
clusters as these evolve with time. Secondly, bright and 
young Hll regions tend to be selected in spectroscopic 
studies of the inner portions of spiral galaxies, while in outer 
disks fainter, and possibly older, nebulae will be typically 
included, because of the r elative scarcity of H ll regions. 
The study of NGC 4625 by [ Goddard iTaL] dSoTT) , however, 
left open the possibility that outer disk Hll regions have 
systematically lower excitation levels relative to inner disk 
ones, which could imply systematic differences in the upper 
initial mass functions. In order to shed light on this issue, we 
have compared our excitation diagrams with the predictions 
from theoretical models of evolving H II regions. 
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Fig. 5. — The [Nii]/Hq: vs. [Oiii]/H/3 diagnostic diagram for the objects shown in Fig.|4](with additional Hll regions from 'Ferguson et al.'1998\ where the 
color code maps the [O III] A50()7/[O II] A3727 line ratio, a measure of the ionization parameter. The model grid by Dopita et al. 1 2UU6 1 tor log(^) = -2 is plotted 



for ages between 0. 1 Myr (thick grey line) and 3.5 Myr, and metallicities Z = 0.2 Zq , 0.4 Zq , 
Outer disk H II regions (R > R25) in different galaxies are represented by white dots. 



Zq and 2 Zq . Dashed lines connect models of constant metallicity. 



4.1. Comparison with model predictions 

The region in the N2 vs. 03 diagnostic diagram that lies 
below the main excitation sequence has seldom received at- 
tention in extragalactic H II regions studies, possibly because 
of the relatively low number of objects found in it. Our inves- 
tigations of outer spiral disks uncovered a significant popula- 
tion of H II regions exhibiting excitation properties that devi- 
ate from the extragalactic H II region excitation sequence. In 
order to clarify the nature of this population we take a look 
at theoretical predictions of the evolution of nebular emis- 
sion line ratios with ionizing cluster age. We consider the 
models presen ted by Dopita et al. (2006), in which the Star- 
burst99 code ( jLeitherer et al.„1999j was used to follow the 
evolution of ionizing clusters syn thesized using the G eneva 
stellar tracks with high mass loss ( Meynet et al.|1994| , and a 
stellar atmosphe re prescription thait i ncludes the non-LTE O- 
star models of Pauldrach et al. (2001) and the Wolf-Ray et star 
models of Hillier & Miller ( ^1998 ). The dynamical evolution 
of the surrounding H II regions was followed by accounting 
for the mechanical energy input from stellar winds and su- 
pernova explosions, and was parameterized in terms of the 
ratio R = (Met / Mq) / (P / k) between the ionizing cluster mass 
and the pressure of the surrounding interstellar medium (the 
latter term in units of 10'*cm"^K). In these models the ion- 
ization parameter, instead of being treated as a free parame- 
ter, is determined by the evolutionary status of the expanding 
H II regions for given values of R and metallicity. The calcu- 
lation of the mo del nebular spectra was carr ied out with the 
Mappings code ( [Sutherland & Dopita||1993] ). For a compar- 



ison with our empirical line flux ratios, we selected models 
with \og(R) = -2, which a re appropriate f or outer disk cluster 



masses of ^10 Mr: 



(Goddard 



et al.|2C 



2010 2011) and an inter- 
We note, however, that 



stellar medium pressure ~ iU" cm 

our conclusions will be largely unaffected by the particular 



choice of R (Fig. 3 of Dopita et al.| 2006 shows that varia- 
tions of R by orders of magnitude affects mostly the ages de- 
rived from the model grids, but not their relative distribution). 
We also remind the reader about the presence of systematic 
uncertainties in the predictions obtained from photoioniza- 
tion models of synthesized stellar population. The calculated 
[O lll]/H/3 ratio appears to be systematically underestimated, 
relative to observations ( [Dopita et al.||2006| attribute the dis- 
crepancy to the inadequacy of the ionizing photon output pre- 
dicted by the stellar atmosphere codes). However, while these 
uncertainties affect the quantitative predictions of the models 
to a certain degree (such as the H II region ages and metallici- 
ties), the qualitative trends are regarded as being quite robust. 
The comparison between the \og(R) = -2 model grid pre- 
dicted by D opita et al. (2006 ) for the N2 vs. 03 diagnostic dia- 
gram at metallicities 'Z = 0.2 Zq, 0.4 Zq, Zq and 2 Zq and our 
extragalactic H II region sample is shown in Fig.lS] In addition 
to the observational data p resented in Fig. | 4| we h ave included 
the H II region sample of Ferguson et al. ( 1998| l, which com- 
prises a handful of nebulae in the outer disks of the galax- 
ies NGC 628, NGC 1058 and NGC 6946. We mapped the 
variation of the [O III] A5007/[O ll] A3727 Une ratio, which re- 
lates directly to the ionization parameter, using the color code 
found at the right hand of the figure. As expected, this ra- 
tio decreases from its maximum value at the top of the dia- 
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gram, along the excitation sequence towards the bottom right. 
We show models for ionizing cluster ages of 0. 1 Myr (thick 
grey line), 1 Myr, 2 Myr, 3 Myr and 3.5 Myr. Models of 
constant metallicity are connected by the dashed lines. The 
match between the model grid and the observational data is 
satisfactory, although objects at the top of the diagram would 
b e better described by models with higher R values (see Fig. 3 
in Dopita et al.|2006| . This comparison indicates that H II re- 
gions older than 3.5 Myr are absent in the data samples we 
included in this diagram. The decrease of the ionization pa- 
rameter with increasing age is well illustrated by the decrease 
of the [O III] A5007/[O ll] A3727 line ratio along tracks of con- 
stant metallicity (the plot also suggests a dependence of the 
ionization parameter on metallicity). 

The interesting point raised by Fig. [5] concerns the region 
of the diagram lying below the main excitation sequence, and 
the distribution of outer disk H II regions in the diagram. The 
models illustrate how, as H II regions of the upper branch of 
the excitation sequence evolve with time, they move to lower 
03 values (at roughly constant N2 values), populating the re- 
gion below the excitation sequence in about 3 Myr. At older 
ages the expanding H II regions have reached such a low sur- 
face brightness that their inclusion in spectroscopic samples 
becomes unlikely. More importantly, H II regions that are lo- 
cated at radii R > R25 in their host galaxies, represented by 
white dots, span the full range of ages from to 3-3.5 Myr, 
with a significant number of them located in the 'zero-age', 
upper portion of the diagram. This indicates that outer disk 
H II regions do not have systematically small [O lll]/H/3 line 
ratios compared to inner disk counterparts. The fact that many 
of them (including a few inner disk regions) are found in the 
lower part of the diagram appears to be consistent with the 
selection effect mentioned earlier, whereby inner disk Hll 
regions are preferentially observed when young and bright, 
while in the outer disks the observer's attention is also drawn 
by older nebulae. 

In order to further test this interpretation we looked for 
trends in the Ha-to-FUV flux ratio, fna/ fpuv, across the N2 
vs. 03 diagram, since this ratio is expected to decrease with 
time, due to the different typical timescales for the (nebular) 
Ha (< 10^ yr) and the (steflai-) FUV (--10^ yr) emission of 
young star clusters. The comparison between Ha and FUV 
emission is not straightforward on a cluster-by-cluster basis, 
because of uncertainties in matching sources at the two dif- 
ferent wavelengths (i.e. due to different pixel scales between 
ground-based optical and spacecraft UV images, but also be- 
cause the FUV emission is intrinsically more diffuse than the 
Ha emission, so that objects that are clearly distinct in Ha can 
often be seen as single clumps in FUV). Therefore, trends in 
the fna / fpuv ratio are expected to be rather crudely defined, 
and only meaningful in a statistical sense. 

For our test we focused on NGC 3621, which hosts a large 
fraction of H II regions falling outside the main sequence in 
the N2 vs. 03 diagram. FUV fluxes were meas ured from 
the available G ALEX images, using SExtractor (Berlin & 
Arnouts 1996| l. We then divided the Hll region sample in 
three broad categories, based o n the age we d erived for each 
object using the model grid by |Dopita et al.|(2006j adopted 
earlier: f < 1.5 Myr, f= 1.5-3 Myr and / > 3 Myr. Since we 
are looking for a statistical trend, and because of the uncer- 
tainties mentioned above, this approach is desirable over one 
with a finer subdivision in age groups. The mean fHa/ fpuv 
values in these three age bins show a decreasing trend with 
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Fig. 6.— The mean Ha-to-FUV flux ratio for the H O regions in NGC 3621 
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age, as illustrated in Fig. [6j as expected. In this figure, we 
normalized the fna/ fpuv ratio to the central bin value, and 
included error bars to indicate the standard errors of the mean 
values. Although the significance of the result is only modest, 
given the size of the error bars, there is a significant difference 
in the fna/ fruv ratio between the youngest and the oldest age 
bin. This ratio decreases by a factor of ^2, in good quantita- 
tive agreem ent with model predictions (see Fig. 7 of Goddard 
20TT) 



et al. 



This result reinforces our conclusion that age 

plays an important role in the distribution of H II regions in 
the N2 vs. 03 diagram below the main sequence. 

An apparent inconsistency in this picture concerns H II re- 
gions hosting Wolf-Rayet (W-R) stars. Since these stars make 
their appearance in the ionizing clusters 2-3 Myr after the ini- 
tial episode of star formation, nebulae containing such stars 
should be located below the main excitation sequence accord- 
ing to the theoretical predictions. However, they are not. As 
an example, the two H II regions in NGC 3621 where we de- 
tected the W-R A4650 feature (numbers 53 and 64 in Table |6]l 
are both located on the excitation sequence. A clue to a pos- 
sible explanation is provided by Fig. |4] (right), where object 
53 is represented by the second largest open circle (i.e. it pos- 
sesses the second largest observed Ha flux), in the upper part 
of the diagram. The presence of W-R stars appears to affect 
the location of the H II regions in this diagram by boosting 
the ionizing flux output of the central clusters, and conse- 
quently providing a larger [O lll]/H/3 ratio. The problem, of 
course, is that this effect should be reproduced by the models, 
which include the W-R evolutionary phase. For a metallic- 
ity Z = 0.4 Zq this phase should last ^2 Myr, which is long 
enough to be resolved by the model grid. 

An important additional point to make concerns the effects 
of metallicity on the evolution on the N2 vs. 03 diagram. At 
high metallicity {^Zq and above) the evolution of the H II re- 
gions occurs along the nearly-vertical, low-excitation arm of 
the main sequence, which is found at [N ll] A6583/Ha ~ -0.4 
(we remind the reader that the [N Il]/Ha ratio is essentially 
a measure of metallicity). Therefore, if the metallicity of an 
aged outer disk H II region approaches the solar value, its lo- 
cation in the diagnostic diagram will remain on, or near, the 
main excitation sequence. This appears to be the case for the 
galaxies M83 and NGC 1512, where very few objects lie out- 
side the region of the diagram enclosed by the curves drawn 
in Fig. |4](we must add that in the case of NGC 1512 our in- 
terpretation is hampered by the small number of outer disk 
H II regions with available [N ll] lines). For a galaxy of lower 
metallicity, such as NGC 3621, we see H II regions straddling 
the space to the left of the near- vertical portion of the main ex- 
citation sequence, starting from [O III] A5007/H/3 ~ 0.6 down 
to ~ -0.4. We conclude that the differences found in the N2 
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vs. 03 diagnostic diagram of outer disk H II regions of dif- 
ferent galaxies stem from differences in metallicity. We also 
infer from the inspection of Fig. [5] that we do not need to in- 
voke a different upper initial massrunction for outer disk pop- 
ulations to explain their distribution in the excitation diagram, 
a possibility we could not rule out in our study of NGC 4625 
( [Goddard et al.|2 01 1 ), but that it is sufficient to account for the 
evolution of the emission line ratios as the ionizing clusters of 
the H II regions become older. 

5. CHEMICAL ABUNDANCES 

5.1. Abundance diagnostics 

For the determination of the oxygen abundances ('metal- 
licities') of the target Hll regions we have used a variety of 
emission line diagnostics. In this way, we expect to be able 
to recover robust qualitative radial trends, regardless of the 
well-known, but unexplained, systematic uncertainties, which 
can amount to a factor of up to ^0.7 dex, presently afflict- 
ing the determination of the chemical abundances of extra- 
galactic H II regions from nebular spectroscopy (see Bresolin 
2008] presoUn et al.||2009al [Kewley & Ellison|[2008| for re- 
cent discussions). In fact, a well-estabUshed, but poorly un- 
derstood, discrepancy exists between the nebular abundances 
obtained from the classical, direct T^-based method ( Menzel 
|et al.|[1941) and the t heoretical predic tions of photoioniza- 



tion model grids (e.g. McGaugh 



1991]l. Various strong-line 



abundance indicators have been introduced to obtain the oxy 
gen abundances of ionized nebulae from the strongest emis- 
sion lines present in their spectra, e.g. when the faint auroral 
lines required in the direct method remain undetected. Cal- 
ibrations of these methods in terms of metallicity have been 
proposed both from empirical abundances (based on auroral 
line detections) and from the results of theoretical photoion- 
ization models. Consequently, the chemical abundances de- 
rived from strong-line methods display large systematic dif- 
ferences when applied to the same observational data, with 
empirically calibrated diagnostics providing measurements at 
the low end of the metallicity range. Differential analysis 
methods are considered to be largely unaffected by these un- 
explained systematic discrepancies, although different cali- 
brations of strong-line methods can yield small differences in 
the slopes of galactic abundance gradients (Bresolin 2011a). 

The faint [O III] A4363 auroral line, which is used, in com- 
bination with the nebular [O III] AA4959, 5007 lines, to ob- 
tain the electron temperature of the ionized gas in the direct 
method, was detected in six Hll regions in NGC 1512 and 
in 12 Hll regions in NGC 3621. These detections are im- 
portant to verify the consistency of the results obtained from 
the remaining strong-line metallicity diagnostics, which do 
not provide direct measurements of the electron temperature. 
Because the strength of [O III] A4363 relative to [O III] A5007 
decreases with metallicity, most of the H II regions in which 
the auroral line was measured are located in the outer disks, 
where metallicities are lower than in the inner disks. The 
direct oxygen abundances were derived from the observed 
line strengths with the nebular package av ailable in IRAF, 
adoptin g the updated atomic data used by Bresolin et al. 
( |2009a| l. We summarize in Table[3]the measured [O III] A4363 
line fluxes and the derived electron temperatures, together 
with the O/H and N/O abundance ratios. The object identi- 
fication numbers in Column (1) refer to those found in the 
catalog of line fluxes presented in Tables [5] and [6] 

In addition, three strong-line nebular abundance diagnostics 
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Fig. 7.— The R23 = ([OlI]A3727 + [Oin]AA4959,5007)/H;3 ratio as a 
function of N2= [Nil] A6583/Ha for the sample of Hll regions shown in 
Fig.H Colored symbols are used for H II regions in the inner and outer disks 
of NGC 1512 and NGC 3621. HlI galaxies extending to low metallicities 
have been added (triangles). The metallicity scale at the far right refers to 
12 + log(0/H) values calculated from the N2 diagnostic. A few outlier Hll 
regions are identified. 

have been used: 

(a) the N2 = [N ll] A6583/Hq; ratio, as calibrated as a function 
of metallicity by Pettini & Pagel ( 2004) from a sample of H II 
regions with [O III] A4363 detections. This diagnostic, there- 
fore, provides oxygen abundances that are approximately on 
the same scale as the auroral line method. 

(b) the R23 = ([O II] A3727 - I- [O ii i] AA4959, 5007)/H/3 ratio, 
introduced by Pagel et al. (1979 1. Of the various R23 cal- 
ibrations as a function of metallicity av ailable in the liter - 
ature, we selected the one proposed by [McGaugh] ( |1991[ ), 
which is based on a grid of photoionization models, and 
which accounts for changes in the nebular ionization param- 
eter through the value of the [Oll]/[Olll] line ratio. It is 
well-known that this calibration, together with other theoreti- 
cal calibrations of R23, provides systematically higher oxygen 
abundances relative to empiri cal calibration s obtained from 
[O III] A4363 detections (e.g. BresoUn et"aL][2009 a), but this 
offset is roughly independent of metallicity, so that relative 
metallicity trends remain well determined. To address the dif- 
ficulty represented by the fact that two different values of O/H 
can correspond to the same R23 value, which gives rise to the 
notion of a 'lower' and an 'upper' branch of the calibration, 
we have plotted in Fig. [7] the R23 indicator as a function of 
N2, which is instead monotonic with O/H. In order to pop- 
ulate the diagram towards the low-metallicity end, we have 
considered, in addition to the H II region samples akeady pre- 
sented in the prev ious figu res, the H II galaxy samples of Te^ 
[levich et aL] ( 1991[ ) and Guse va et"ar] ( |201l| l, shown by the 
triangle symbols. It can be seen that for log(N2) values down 
to at least -1.5, corresponding to 12-i-log(0/H)~8.1 in the 
empirical abundance scale (shown to the right of the H II re- 
gions sequence), we can safely select the upper R23 branch. 
All of the objects we observed in NGC 1512 and NGC 3621 
fall in this range, and we have therefore applied the analyti- 
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cal expression for t he upper branch calibration of |McGaugh 
(1991 ) provided by lKuzio de Naray et aLl ( |2004l ). ITie plot in 
Fig. n\ identifies a few outliers relative to the H II region se- 
quence. To the right of the sequence, in particular, we find 
nebulae with particularly hard ionizing spectra, which can be 
related to the presence of He II A4686 emission (e.g. region 1 
in NGC 3621; for other regions fa lling in this category and 
not shown here see |Bresolin |2011bi), or shocked gas si gna- 
tures (e.g. NGC 4395-001 -(J05, see |McCall et al.|1985| l. 

(c) th e empirical method proposed by Pilyugin et al. 
(20101, which calibrates the strengths of the [O llj and 
[O III] lines as a function of the O/H ratio, based on 
[O III] A4363-based abundances, but making use also of the 
[N II] AA6548, 6583/H/3 and [Sll] AA6717,6731/H/3 line ra- 
tios to empirically constrain the electron temperature. This 
method is being referred to as the ONS method. Once again, 
by construction this diagnostic should provide metallicities 
that are consistent with those derived from the [O III] A4363 
line. Compared to previous empirical calibrations provided 
by the same authors, the ONS method has the advantage of 
being applicable to the full range of metallicities encountered 
in extragalactic H II regions. 

5.2. Radial metallicity gradients 

The galactocentric radial gradients of the O/H abundance 
ratio are shown in Fig. [8] for NGC 1512 and Fig. [9] for 
NGC 3621. The upper panels present the abundances ob- 
tained from the [O III] A4363 diagnostic (red-yellow circles) 
and N2 (open circles). The lower panels show the abundances 
obtained from R23 (open triangles) and the ONS method (open 
circles). Linear regression fits to the inner disk data {R < R25) 
are represented by the sloping grey lines. The different diag- 
nostics show that the oxygen abundance distributions are flat 
beyond /?25- Within the errors, gradient determinations in the 
outer disks are consistent with a flat radial trend. For sim- 
plicity, we have then represented the outer disk metallicities 
with the average O/H ratios, shown by the horizontal lines. 
The statistical test described below indicates that the change 
in slope takes place near the isophotal radius, in the radial 
range 0.8 - 1.2 /?25- We did not attempt to pinpoint the exact 
radius at which this occurs, since a significantly higher num- 
ber of targets in the transition zone would be required. The 
small discontinuities in O/H seen in Fig. [8] and l9] reflect the 
difficulty to establish the radius at which the transition occurs. 

The statistical significance of the flattening of the abun- 
dance gradient beyond the isop hotal radius has been e stimated 
with an F-test, following Za hid & Bresolin] ( |2011[ ) in their 
study of NGC 3359. The test allows us to check whether 
the improved fit (i.e. smaller residuals) obtained using a two- 
slope model for the radial abundance distribution (in which 
we set the outer gradient to be flat, for the reason given above) 
is statistically significant compared to a single-slope represen- 
tation. In general we find that the flattening is significant at 
the >99.99% level. The exception is represented by the gra- 
dient observed with the N2 diagnostic in NGC 1512, where 
a single slope (dashed line in Fig. IHll reproduces the observed 
data as well as if we used two lines of different slopes. How- 
ever, this is a result of the small number of H II regions in 
the outer disk of this galaxy for which we can calculate N2- 
based abundances (see below). The R23 metallicity indicator, 
for which we have more data points in this galaxy, displays a 
statistically significant flattening. 

In the following we discuss separately our results for 



NGC 1512 and NGC 3621. 

NGC 1512 

As mentioned in Sect. |2.2| we do not have red spectra 
covering the Ha and [Nil] lines for the 10 Hll regions lo- 
cated in the outer eastern spiral arm, and therefore the N2 
and ONS diagnostics provide an incomplete picture of the 
abundance distribution in the outer disk of NGC 1512. How- 
ever, the existing data are in good qualitative agreement with 
the R23 abundances, which clearly show a flat metallicity be- 
yond the isophotal radius. This is also confirmed by the few 
[O III] A4363 detections, even though the O/H value for the 
outermost data point, corresponding to region 61, is slightly 
below the rest (although not significantly, considering the er- 
ror bars; see upper panel of Fig. [8]l. The N2 and the auroral- 
line based abundances are in good quantitative agreement for 
objects in which both methods could be used. 

The abundances obtained from the N2 diagnostic reveal 
that five inner disk H II regions, two of which have consis- 
tent [O III] A4363-based O/H abundances, lie significantly be- 
low the regression line. Interestingly, all of these objects, en- 
closed in the box drawn in Fig.|8|(top) and numbered 53 to 
57 in F ig. \i\ are located in the bridge' structure identified 
by [Kori balski & Lopez-Sanchez (2009) in the zone between 
the centers of NGC 1512 and NGC 1510, and possibly corre- 
sponding to tidal debris (see Fig. [Til. These H II regions there- 
fore do not appear to belong to the same population of H II re- 
gions found in the main body of NGC 1512, and might instead 
be related to material that was stripped during the interaction 
between these two galaxies. The lower O/H abundance ra- 
tio measured in the bridge is consistent with the idea that this 
gas has reached its current location from larger radii, possibly 
around the isophotal radius (or beyond), where a similar O/H 
value is found. Although this association to tidally stripped 
material remains speculative, we have removed objects 53-57 
from the calculation of the linear regression of the metallic- 
ity as a function of radius, parameterized as 12 H- log(0/H) = a 
+ b (R/R25). The results obtained for coefficients of this fit 
from the different abundance diagnostics are summarized in 
Table g] It can be seen that the N2 and the ONS meth- 
ods provide results that are consistent with each other, for 
both the intercept value, 12H-log(0/H)~8.8, and the slope 
of the abundance gradient in the inner disk, approximately 
-0.034 dexkpc"'. The R23 indicator yields a ^0.4 dex larger 
intercept value, as expected, but also a steeper slope, although 
in this case the linear fit displays an rms scatter that is almost 
double the one measured using the N2 diagnostic. We finally 
note that, after discarding the objects in the bridge, the rms 
abundance scatter obtained for the inner disk of NGC 1512 
from the N2 diagnostic, 0.06 dex, is the same value that we 
measure for NGC 3621. 

NGC 3621 

Thanks to the large number of H II region abundances avail- 
able in both the inner and the outer disk of this galaxy, 
the manifestation of a break in the abundance gradient of 
NGC 3621, and the flat distribution beyond the isophotal ra- 
dius, are very well established by any of the strong-line metal- 
licity diagnostics we adopted. The N2 abundances are in good 
agreement with those derived from the auroral line method, in 
those cases where the [O III] A4363 line was measured. As 
in the case of NGC 1512, the N2 and ONS diagnostics yield 
consistent abundance gradients, which are well anchored to 
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Fig. 8. — (Top) The galactocentric O/H abundance gradient in NGC 1512, 
as determined from the N2 diagnostic. Abundances obtained from the 
[O III] A4363 auroral line are shown as yellow dots. The dotted box encloses 
the 5 Hll regions that are located in the 'bridge' between NGC 1510 and 
NGC 1512. The galactocentric distance axis is labeled both in terms of S/R25 
(bottom scale) and in kpc (top scale). The dashed line shows a single-slope 
linear fit for all N2-based data points. (Bottom) The galactocentric O/H abun- 
dance gradient in NGC 1512, as determined from the R23 indicator (open 
triangles) and from the ONS method (open circles). In both panels the full 
grey lines represent the linear fits to the inner disk H II region radial abun- 
dance distribution (R<R2s), and the constant values assigned to the outer 
disk from the mean O/H ratio for R > R25 ■ For clarity, error bars are shown 
only in the top panel. 
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Fig. 9. — Same as Fig.lSlfor NGC 3621. A few outliers discussed in the 
text are labeled by their number in TablelS] 



the empirical [O ill] A4363-based metallicities. The overall 
metallicity is about 0.1 dex lower than in NGC 1512. While 
the slope derived from R23 is again slightly steeper, it is not 
significantly different from the one derived from N2. For all 
the diagnostics, the rms abundance scatter relative to the lin- 
ear fit is small, ^0.06 dex. The abundance distribution is thus 
quite tightly defined, and its scatter is essentially determined 
by the observational errors, rather than true metallicity varia- 
tions at a given galactocentric radius. 

In Fig. |9] we identify a few outliers relative to the radial 
trend defined by the majority of the H H regions. In the top 
panel, region 46 has an N2 -based abundance that is about 
0.2 dex lower than the [O Hi] A4363-based one. The latter 
matches the abundance from the linear fit at the correspond- 
ing radial distance, so apparently there is something peculiar 
about its N2 line ratio. In the bottom panel, the metallicity of 
region 1 estimated from R23 is clearly underestimated , a re- 
sult that is due to its high excitation level. Bresolin^^2011b| l 
showed that this is typical for H 11 regions displaying nebular 
He n emission. There are also three objects (54, 57 and 73) 
for which the ONS method appears to underestimate the O/H 
abundance ratio (these were excluded from the calculation of 
the mean outer disk abundance given in Table Hll. Pi lyugm] 

tal. (2010) provide three separate analytical expressions to 
rive O/H, for different ranges of the [N 11] AA6548, 6583/H^ 
and [S II] AA6717, 6731/H^ line ratios. Owing to the observed 
log([Nll]/[Sll])<-0.25, these Hn regions are classified as 
'hot', and are assigned a rather low metallicity. If they were 
classified as 'warm', by a change in the criterion used to sep- 
arate H n regions in different temperature regimes, the O/H 
ratio for these three H n regions would agree with the rest of 
the outer disk sample. This sugge sts tha t some adjustment 
to the criteria adopted by Pilyugin et al. (2010 1 could lead 
to a better estimation of the O/H ratio at metallicities around 
12H-log(0/H) = 8.2. 

5.3. Nitrogen 



In Fig. 10 we show the O/H vs. N/O diagram for Hll re- 
gions in the outer disks of NGC 1512 and NGC 3621, but 
only for those targets for which we have a firm detection of 
[O 111] A4363, and therefore of the electron temperature (those 
inclu ded in Table [3|). We inc lude targets in th e outer disks of 
M83 dBresohn et al.||2009bD and NGC 4625 dGoddard et al] 
201 l|l. For comparison, we also show a sample of extra^ 



galactic H n regions (triangles;, Garnett et al. 1997; KennicuH] 
et al.-2003a;"Bresolin et al.|20(J5||2009a||Esteban et al |2009r 
Bresolin et al. 2010) and H n galaxies (empty squares; |Izotov| 
et al. 1994; Guseva et al. 2009). In all cases, we only consider 
abundances obtained from the knowledge of the electron tem- 
perature, as derived from the detection of auroral lines. 

The distribution of the points representing the outer disk 
H II regions is roughly consistent with the well-known gen- 
eral trend in this diagram, in which a flat distribution at low 
O/H ratios with constant N/O ~ -1.5, attributed to predom- 
inant primary nitrogen production at low metallicity, is fol- 
lowed by an increase of N/O with O/H, as a signature of an 
additional, secondary production of nitrogen at higher metal- 
licity. Considering the large spread in N/O at constant O/H, 
the fact that outer disk Hll regions on average do not dif- 
fer significantly in a systematic way from disk H 11 regions 
in terms of their N/O ratio excludes gross deviations of the 
initial mass function (IMF) of the embedded star clusters in 
the extended disks from the canonical one, typical of inner 
disk clusters. A hypothetical systematic decrease in the up- 
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TABLE 3 
Direct method electron temperatures and abundances 



ID 


[Olll] 
4363 


(K) 


12 + log(0/H) 


log(N/0) 


(1) 


(2) 


(3) 


(4) 


(5) 



NGC 1512 

1 5.9 ±0.7 12090 ±750 8.21 ±0.11 

12 5.8 ±0.3 11980 ±460 8.18 ±0.07 

45 3.4 ±0.4 11390 ±610 8.17 ±0.10 

55 2.3 ±0.2 11440 ±540 8.16 ±0.08 

57 2.4 ±0.2 10810 ±500 8.23 ± 0.09 

61 6.9 ±0.6 13380 ±710 8.05 ± 0.09 

NGC 3621 

1 15.2 ±1.2 13980 ±740 8.21 ± 0.08 

16 1.7 ±0.2 10360 ±550 8.25 ± 0.09 

20 1.7 ±0.2 10060 ±450 8.30 ± 0.08 

46 2.0 ±0.1 9250 ±270 8.44 ± 0.06 

49 3.3 ±0.3 11350 ±480 8.19 ±0.08 

51 2.5 ±0.2 10280 ±360 8.25 ± 0.07 

52 1.8 ±0.2 9700 ±400 8.36 ± 0.08 

53 2.8 ±0.1 9890 ±270 8.37 ± 0.06 

60 3.1 ±0.3 10570 ±510 8.33 ± 0.09 

62 2.1 ±0.3 10540 ±570 8.27 ± 0.09 

63 3.5 ±0.2 11250 ±440 8.24 ± 0.07 

64 3.5 ±0.2 11510 ±440 8.18 ± 0.07 

Note. — The reddening-corrected intensity of [O III] A4363 in Column (2) is in units of H/3 = 100. 



-1.32 ±0.18 
-1.35 ±0.13 

-1.35 ±0.12 
-1.26 ±0.14 



-1.31 ±0.14 
-1.47 ±0.13 
-1.18 ±0.12 
-1.30 ±0.09 
-1.50 ±0.11 
-1.34 ±0.10 
-1.56 ±0.11 
-1.48 ± 0.08 
-1.40 ±0.12 
-1.42 ±0.13 
-1.50 ±0.10 
-1.47 ±0.10 



TABLE 4 
Results of linear regression to the abundance gradients 



Diagnostic 


Intercept (a) 


Slope (b) 
dexkpc-' dexRjs"' 


rms 


Outer disk 
mean 


(1) 


(2) 


(3) (4) 


(5) 


(6) 







NGC 1512 








N2 

ONS 

R23 


8.83 ±0.03 
8.75 ±0.03 
9.22 ±0.05 


-0.034 ± 0.004 
-0.033 ±0.005 
-0.059 ±0.007 

NGC 3621 


-0.40 ±0.04 
-0.39 ±0.05 
-0.69 ±0.08 


0.06 
0.07 
0.11 


8.31±0.11 
8.32±0.04 
8.54±0.08 


N2 

ONS 

R23 


8.73 ±0.03 
8.68 ±0.03 
9.09 ±0.04 


-0.054 ±0.006 
-0.047 ±0.006 
-0.066 ±0.007 


-0.51 ±0.05 
-0.44 ±0.05 
-0.62 ±0.06 


0.06 
0.06 
0.07 


8.23±0.05 
8.23±0.05 
8.59±0.08 



per mass cutoff would lead to an increased N/O ratio, due 
to the removal of the main producers of oxygen, i.e. massive 
hot stars. While we postpone a more detailed discussion to 
a future work (Bresolin et al., in preparation), we exclude a 
dramatic reduction of the upper mass cutoff from 100 Mq to 
25-30 M0. 
However, Fig. 



10 suggests that outer disk Hll regions in 



NGC 3621 lie towards systematically low N/O values, com- 
pared to similar objects we studied in M83, NGC 4625 and, 
in part, NGC 1512. While this might be simply due to small 
number statistics, it is tantalizing to speculate in terms of the 
different evolutionary status of the extended disks of these 
galaxies. Among the four galaxies, NGC 3621 is the only 
one which we can consider isolated and non-barred. The oth- 
ers show clear signs of tidal interactions (as in the case of 
NGC 1512 and NGC 4625) or mi ght ha ve had an interaction 
in the past (as suggested for M83, van den Bergh 1980), and 
contain stellar bars (NGC 1512, M83). Both of these proper- 
ties are linked to the presence of radial gas flows, which can 



presumably mix metals out of the main disk of the galaxy into 
the outer disk. We can speculate that the extended disk of 
NGC 3621, in the absence of such a radial mixing process, is 
less chemically processed than the remaining extended disks. 
Because of the longer timescale for nitrogen production com- 
pared to that of oxygen, the extended disk of NGC 3621 
would therefore display a lower mean N/O ratio. Of course, 
we realize that observations of the N/O ratio in a significantly 
larger sample of extended disk galaxies will be necessary be- 
fore this hypothesis can be seriously tested. 

6. DISCUSSION 

In the previous section we illustrated the flattening of the 
Hll region radial abundance gradient occurring in the ex- 
tended disks of the spiral galaxies NGC 1512 and NGC 3621. 
The O/H ratio is essentially constant for galactocentric dis- 
tances larger than ^10 kpc, to the outermost regions we 
observed at distances of 18 kpc (NGC 3621) and 25 kpc 
(NGC 1512). The small abundance scatter that we mesaure. 
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Fig. 10.— 0/H vs. N/O for extended disk Hll regions in NGC 1512, 
NGC 3621, M83 and NGC 4625, compared with a sample of extragalactic 
H II regions (triangles) and H II galaxies (open squares) taken from the liter- 
ature (see text for references). Only objects with abundances obtained from 
the knowledge of the electron temperatures (derived from auroral lines) are 
displayed. 



on the order of 0.06 dex, is consistent with what is expected 
from the observational uncertainties. Strong evidence for flat 
outer disk gradient s in external spiral ga laxies has been re- 
cently presented by Bresolin et al. (|2009b|l for the prototypical 



XUV disk galaxy M83, an d subseq uently by Goddard et al. 
(|20n) for NGC 4625, and |Werk e t al. (201jjl for a sample 
of 13, mostly merging or interacting, galaxies, although in 
the latter case the abundance scatter is in general consider- 
ably larger (^0.15 dex, compared to 0.06 dex), due to the 
significantly smaller number of H 11 regions observed. A flat- 
tening, or possibly an upturn, in the metallicity gradient has 
also been detected from stellar photometry of the old, resolved 
stellar populatio n (red gi ant branch stars) in the outer d isks of 
NGC 300 (Vla jic et al.| [2009) and NGC 7793 (Vlaji c et al. ' 



201 1 1, even though this result is somewhat affected by the 



age-metallicity degeneracy. Strong indications for a flat outer 
abundance gradient have also b een found for the Milky Way 
from a variety of studies (e.g. |PedicelU et al.|[2009t |Lepine| 
|et al. [20111 and references therein). A schematic representa- 
tion of the galactocentric O/H abundance gradients obtained 
by our group for the extended disk galaxies M83, NGC 4625, 
NGC 1512 and NGC 3621, determined from the N2 metalHc- 



ity indicator, is presented in Fig. 1 1 This provides an impres- 
sion of the similarity in the behavior of the metallicity gradi- 
ent observed at large galactocentric distances, despite the vari- 
ety in this (admittedly small) sample of galaxies (for example 
in terms of physical size, interaction activity and presence of 
bars). The value of 12H-log(0/H) in the outer disks ranges 
between 8.2 and 8.4, in a way that appears correlated with the 
mean metallicity of the inner disk (as discussed below). 

Negative radial abundance gradients are a weU-estabhshed 
feature of the inner disks of spiral g alaxies ( jVila-Costas &| 
Edmunds 19921 Zaritsky et al.| 1994| i. They are understood 
within the inside-out galaxy formation scenario as a result of 
the increased timescales of the gas infall with radius within a 
A cold dark matter context and the consequent radial depen- 
dence of the star formation rate. On the other hand, observa- 
tional and theoretical investigations of the metal content in the 
outer regions of galaxies have only recently begun. With the 
discovery of UV emission from star-forming regions located 
in the extended outer disks of a significant fraction of galax- 
ies in the Local Universe (Gil de Paz et al. 2005 , Thilker et al.. 
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Fig. 11. — Schematic representation of the galactocentric O/H abun- 
dance gradients obtained by our group for the extended disk galaxies M83, 
NGC 4625. NGC 1512 and NGC 3621. using the N2 metallicity indicator. 
Dots are shown at R = 0.4 Ris for each galaxy to represent the characteristic 
abundances of their inner disks. 



|2005| |2007| l, the crucial role played by studies of the outer 
fringes of galaxies in our understanding of galaxy formation 
and evolution has become evident. While recent investiga- 
tions have focussed especially on the star forming properties 
of these XUV disks (|Dong et al. 2008; Bigiel et al. 2010b| 
[Goddard et al.|2010[ [Alberts et al..2011) , it is essential to in- 
clude the additional constraints to galactic evolution studies 
provided by the chemical abundance analysis, of H II regions 
in particular, although this is hindered by observational diffi- 
culties, due to the relative scarcity and faintness of the targets. 

6.1. The flatness issue 

For the interpretation of our observational results it is im- 
portant to first address the possible reasons for the flat chem- 
ical abundance distribution in extended galaxy disks. A po- 
tential clue is offered by the comparison of the star formation 
efficiency (SFE) be tween inner and outer disks of galaxies. 
[Bigiel et aL] ( |2010b) l showed that this efficiency, defined as the 
ratio of the surface density of star formation rate to that of the 
available neutral hydrogen gas, SFE = J^sFR/^Hh is consider- 
ably flatter in the outer disks, in the radial range between one 
and two isophotal radii, when compared to the inner (R < R25) 
disks. Examples where this behavior continues t o even larger 
galacto centric distance s can be found in M 83 (Bigiel et al. 
[2010a| l and NGC 864 ( [Espada et aLllSoTT) . It is interesting 
to note that the flattening in M83 occurs around R/R25 = 1.5 
(see Fig. 4 of Bigiel et al. "2010a), which c orresponds to the 
approximate radius at which ,Bresolin et al. (2009b) found the 
radial abundance gradient to flatten out (for 1 < R/R25 < 1 .4 
there seems to be a 'superposition' of inner and outer disk 
abundance gradients). 

The gaseous oxygen-to-hydrogen abundance ratio per unit 
disk surface area can be approximated by 



O 
H 






cxSFE 



(1) 



where yo is the net oxygen yield by mass, /i = 11.81 is the 
conversion factor from number ratio to mass fraction, and t is 
the timescale for the star formation activity. This expression 
suggests that a slow radial dependence of the star formation 
efficiency at large galactocentric distances could translate into 
a flattening of the metallicity gradient. 
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A potential complication to this picture is represented by 
the relatively l arge O/H value measured in the outer disks. As 
pointed out by'Bresolin et al.'('2009 b|) i n the c ase of M83, and 
later confirmed by Werk et al. (2010) [201 l| l, outer disks are 
overabundant for their (extremely large) gas fractions. Virtu- 
ally all of the outer disks in which nebular abundances have 
been derived todate possess very extended neutral hydrogen 
envelopes, which act as large gas reservoirs for the ongoing 
low-level star forming activity we witness in the form of UV- 
emitting knots and H II regions. Due to the preponderance of 
H I in the calculation of the baryonic mass in the outer disks, 
the effective yield jgff = Z/(/n/x~') (where fi is the gas frac- 
tion) is thus quite high. 

Considering the O/H values anchored to the [O III] A4363 
auroral line-based abundances as lower limits to the nebu- 
lar oxygen abundance (since the various strong-line methods 
provide equal or higher O/H abundances), we see from Ta- 
blefflthat 12H-log(0/H) is > 8.2 and > 8.3 in NGC 3621 and 
NGC 1512, respectively. An estimation of the time necessary 
to enrich the outer disks to these abundance levels, assuming 
a constant star formation rate (SFR) equal to the presently 
observed value, can be quickly obtained from the knowl- 
edge of the SFR and the surface density of gas (see Eq. 1). 
While in this way we obviously neglect the (unknown) de- 
tails of the chemical enrichment history of the galactic disks, 
the approximation will provide an idea of the timescales in- 
volved. The calculation is done first for the case of the outer 
disk of NGC 3621, for which the relevant rn ean quantities 
are readily available from Bigiel et al. (2010b): T,sfr = lO"'* 
Mq yr~' kpc~^, S^/ = 5 MqPC"^. Adopting an oxygen yield 
of 0.01 ( Maeder|1992| l, the time required to enrich the inter- 
stellar medium up to 12H-log(0/H) = 8.23 is --10 Gyr. The 
timescale would more than double if we wished to obtain 
12 H- log(0/H) = 8.6, as measured from the R23 abundance in- 
dicator. 

Even though a timescale that is consistent with the age ex- 
pected for inner galaxy disks (> 10 Gyr) has been found in 
this casa^ outer disks are unlikely to have formed stars for 
such long times. In the inside-out scenario for galaxy growth 
the outer disks of galaxies are considerably younger. Cos- 
mological hydrodynamical simulations indicate that the outer 
portions of galaxy disks could have mean formation times 
around 4-6 G yr ago, while the inner few kpc are formed > 10 
Gyr ago (e.g. Scannapieco et al.|2008 2009 1. Supporting ob- 
servational evidence has recently been provided by the star 
formation histories of nearby resol ved spiral galaxies, such as 
M33 (Williams et al.|2009;.Barker efal.pOl l| l. IGogart en et al.| 
( |2010p found that while for the inner 3 kpc of NGC 300 50% 
of the stars had already formed 1 1 Gyr ago, at a galactocentric 
distance of 5 kpc (corresponding approximately to the isopho- 
tal radius) the same level of star formation was reached about 
6 Gyr ago. Age gradients (younger stars at larger radii) have 
also been inferred from resolved stellar populations studies 
in the outer disks of spiral galaxies such a s M81 (Davidge] 
|2009] l and the extended disk galaxy M83 ( JDavidgeipcnnT 
ITie absence of a significant populati on of r ed giant branch 
stars in the outer disk of M83 noted by Davidge (2010) points 
against a significant star formation activity a few Gyr before 
the present time. Hl-rich galaxies with star-forming outer 



' Our calculation for M83 in [Bresolin et al.|J2009b| contained an error, 
which led to underestimate the timescale by a tactor ot jj. = 11.81. Therefore, 
for the outer disk of M83 the required timescale for chemical enrichment 
would be ~3 X lO'" yr, much longer than a Hubble time. 



disks have been shown to have color gradients that are consis- 
tent with rec ent gas accretion and an inside -out growth (Wang] 
et al.[[20lT i. The sub-critical density of the gas available in 
the outer disks for star formation suggests that it is unlikely 
that significantly higher star formation rates have been expe- 
rienced in the past. We therefore argue that the systems we 
are analyzing are observed in a phase that could have started 
a few Gyr ago (consistent with the age estimates obtained for 
the oldest UV-bright clusters), during which a large amount 
of neutral gas has become available for star formation, and 
during which the mean star formation rate has remained close 
to the present-d ay value (as als o suggested by the outer disk 
simulations by Bush et al.|2010) . The problem of the chemi- 
cal enrichment of outer disks is exacerbated under these con- 
ditions, because of the reduced amount of time available for 
metal production. For example, in NGC 3621, star forma- 
tion at the currently observed levels would enrich the ISM to 
12H-log(0/H)~7.8in4Gyr 
For the case of the outer disk of NGC 1512 we have de- 



rived the following appro ximate parameters from Koribalski 
& Lopez-Sanchez ([2009[): Esfr = 2.5 x 10"^ M© yr"' kpc"^ 



Uh/ = 2 M0PC . ITie resulting timescale to reach the O/H 
abundances in Table[4[is ^2-3 Gyr, depending on the metallic- 
ity diagnostic used. The ongoing interaction with NGC 1510 
might explain the shorter value for the estimated timescale 
compared to NGC 3621 (the star formation rate could have 
been recently enhanced, thus leading to overestimate the aver- 
age SFR over time), although this remains rather speculative. 
We conclude that, at least for a few galaxies we have 
analyzed so far (such as M83 and NGC 3621), it is un- 
likely that the in situ star formation could have produced and 
ejected enough metals to enrich the interstellar medium to 
the presently observed values. We also note that O/H self- 
enrichment of the H II regions from embedded massive (O 
and Wolf-Rayet) stars and type 11 supernovae appears a very 
unlikely mechanism, especially at metallicities below solar, 



in ligh t of the theoretical results presented by Kroger et al. 
( [2006| and |Wofford| ( [2009] l. 

6.2. Radial metal mixing 

Although we should expect a fraction of the elevated metal 
amounts found in the outer disks to be produced in situ, since 
stars are known to populate galaxy disks to large galactocen- 



tric distances ( [Bland-Hawthorn et al. 2005 [ Davidge 20101, 
another considerable fraction must have reached the outer 
disks, possibly from the inner disks, via some transport mech- 
anism, efficient out to galactocentric distances of at least 20 - 
25 kpc, a s observed in NGC 15 12, NGC 3621 (this work) 
and M83 (Bresolin et al.ll2009bl). Several possibilities have 
already been poin ted out by other authors ( Werk et al.|2010 
20TT| |Vlajic"eFal.|2009| |201 Ip . These include the effects~' 



radial gas flows induced by bars ([Fri edli et al. 11994) and tidal 
interactions, as already discussed in Bresolin et al. (2009b). 
Abundance gradient breaks have been reported in a few barred 



galaxies (e.g. NGC 3359, Martin & Roy 1995 Zahid & 



BresoUn|201 1 1, although it should be noted that such breaks 



are expected to occur near the bar corotation radius, in most 
cases located within the galaxy optical disk, so that this would 
not in general explain breaks taking place further out around 
the isophotal radius. The flattening of radial abundance gradi- 



ents d ue to galaxy tidal interactions and merging (Rupke et al. 
I could explain most of the cases investigated by Wer^ 
( [201 1| , as wefl as our result for NGC 1512. However, 



2010 
eTaT 



the fact that apparently isolated galaxies (although we cannot 
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exclude past merger events), such as NGC 3621, also display 
similar uniform metallicity distribution in their outer disks 
leads us to conclude that alternative mechanisms of metal re- 
di stribution must play an important role. 

|Werk et~aL |20T1) argued in favor of the mixing sce- 
nario without feedback and outflows proposed by Tassis et aL] 
(12008 ), in which metals are driven in the hot gas phase to large 
galactocentric distances by turbulence and large-scale mo- 
tions resulting from gravitational instabilities, accretion flows 
and mergers. This model was developed to ex plain known 
scaling relations of dwarf galaxies. According to Tassis et aT] 
(J2008 ), mergers should be the dominant source of turbulence, 
tnerefore it is unclear whether this mechanism can fully ex- 
plain the chemical enrichment of the outer disks of nearby, 
non-dwarf isolated galaxies. The model seems also incompat- 
ible wi th the direct evid ence for gas outflows in dwarf galaxies 
(e.g. Elson et al.|2"0TT] l. 

Additional processes of angular momentum transport that 
can help to explain the ra dial mixing of metals in galaxy disks 
include radial gas flows (Lacey & Fall'1985', 'Goetz & KoepJ 
pen 1992; Portinari & Chiosi 2000; Ferguson & Clarke 200Tf 
Schonrich & Binney 2009 ; Spitoni & Matte u cci|201 1[ among 
others), perturbations by satellite galaxy encounters ( Quillen| 
leTal. 2009; Bird et al. 2012), reso nance scattering w ith tran- 
sient spiral density waves ( ,Sellwood & Binneyl20 02[l, and th e 
overlap of spiral and bar resonances (Minchev e t al.||20l"T| l. 
Stellar radial migrations ("Roskar et al. 2008b) can lead to a 
flattening of the metallicity gradients for old stellar popula- 
tions with increasing radius (Roskar et al. 2008a), as observed 
for example in NGC 300 ( Vlajic et al._2009) . However, since 
the effect of radial mixing on age and rnetallicity as a function 
of radius is cumulative with time, and since stars and gas are 
essentially decoupled, we do not expect the same mechanism 
to be able to explain the flat abundance gradients we derive 
for the gas out to galactocentric distances of ^20 kpc. 

6.3. Enriched infall 

An alternative interpretation of the chemical properties of 
extended disks involves the predictions from recent cosmo- 
logical simulations in the hierarchical context concerning the 
evolution of galaxies that consider cold gas accretion, mod- 
ulated by feedback , outflows and star forma t ion processes 
(|Keres et al.|[2005) IDekel & Bimboiml [20061 [Ocvirk et al.| 



2008 Bouche et a 



2010|). The direct detection of gas ac- 



cretion into galaxies is elusive, and has been inferred from 
gas depletion timescales, which are considerably shor ter than 
the H ubble time in galaxies at low and high redshift ( Genzel| 
|et al.p OlO), and from considerations concerning misaligned 
components in disks, such as warps, polar disks and other 
ring structures (Roskar et al. 2010 Finkelman et al. 2011 
Spavone et al.|2011V Metal-poor Lyman limit systems have 



been recently proposed as direct observational support for 
cold stream accretion (IRibaudo et al. 11201 11 IFumagalU et al.l 

|2gTT) . — 

G alactic scale outflows are well established observationally 
(e.g. Veilleux et al.||20(J5 Steidel et al.||2010 and references 
therein). The ejection of chemically enriched gas from galax- 
ies appears necessary to explain the enrichment of the circum- 
galactic an d intergalactic medium via galactic winds (Oppen^ 
[heimer &D ave 2008 ; Tumlinson et al. 201 1 ), the origin of the 
mass-metallicity relation (Finlator & Dave|2008| ), and to cir- 
cumvent discrepancies between predict ed and observed b ary- 
onic mass fractions locked in galaxies ("Kere s et al.|2009 1. In 
recent outflow models a 'wind recycling' process ensues, in 



which the infalling gas is progressively enriched by metals 
produced by stars within galactic disks, ejected by outflows 
over distances of ~80 kpc and reincorporated into the parent 
star-forming galaxies over tim escales on the order of 1 Gyr 
( [Oppenheimer & Dave 2008). This process was more effi- 
cient at early epochs and high redshift, z > 2, when accretion 
and the star formation rate were significantly higher than at 
z = 0, giving rise to rapid metal enrichment in the Universe. 
In the present epoch the chemically enriched material located 
in the intergalactic medium participates in ongoing, weaker 
gas inflows. Within this picture, even galaxies that are not 
massive enough to sustain important outflows can thus accrete 
enriched gas. 

The metallicity expected for infalling gas in the local Uni- 
verse does not necessarily have to be 'primordial' or very 
low (0. 1 Zq or below). High velocity clouds represent an 
important mode of gas accretion for the Milky Way, and are 
also observed in connection with other g alaxies (e .g. ISancisi 
[et al. 2008 ). Cosmological simulations by Keres & Hemquist 
^009 j ) show how such clouds can be produced from the fila- 
mentary structure in cold gas accretion, providing Milky Way- 
size galaxies with a continuous source of gas for star forma- 
tion. Absorption line measurements indicate metallicities in 
the range 0. 1 Zq to supersolar ( Wakker 2001 ; Yao et al. 201 1 ). 
This range probably reflects the variety in the origin of these 
systems, from metal-poor intergalactic medium accretion and 
small me rgers to met al-rich gas participating in the galactic 
fountain ( [Shapiro &~F ield 1976). Whatever their origin, en- 
riched gas clouds (0.2 - 0.5 Z©) associated with the Milky 
Way can be found at galactocentric di stances (9-18 kpc) t hat 
place them in the outer Galactic disk (Tripp & Song|20l"T) 



The recent cosmological hydrodynamic simulations by 



Dave et al. (201 lb I highlight the importance of reaccretion 
of metal-enriched gas in defining the observed galaxy mass- 
metallicity and mass-gas fraction relations at redshifts z < 1 . 
These authors followed the redshift evolution of the gas- 
phase metallicity (in terms of the star formation-weighted 
O/H ratio) of both the interstellar medium of model galax- 
ies and the infalling gas. Their momentum-conserving wind 
models at z = reach an ISM metallicity approximately be- 
tween 12H-log(0/H) = 8.3 (2 X lO'^M© galaxy mass model) 
and 12H-log(0/H) = 8.6 (10'°Mq). The corresponding infall 
gas metallicities are roughly 0.15 dex lower, with the ratio 
between the infalling gas metallicity and the inner disk ISM 
metallicity az — 0.7. This parameter enters directly into the 
equation for the equilibrium metallicity of the gas under the 
combined eff ects of star formation and dilution by infall (see 
[Dave et al.|2 011b a for details). If we consider the abundance 
gradient plots for NGC 1512 (Fig.[8]l and NGC 3621 {¥igM 
we see that the difference between the mean O/H ratio for the 
inner disks and the outer disk is comparable to this theoret- 
ical value. If we use the inner disk metallicity measured at 
/? = 0.4/?25 as representative of the 'effective' or characteris- 
tic metallicity ( [Zaritsky et al. 11994] see the dots in Fig.fTT]) we 
obtain, from the results for the N2 metallicity indicator in Ta- 
ble[4| az — 0.5. Interestingly, if we also consider the extended 
diskgalaxies M83 and NGC 4625 we studied previously, we 
obtain similar values {a = 0.46 and 0.43, respectively). In this 
limited sample of extended disk galaxies we thus find that 
the outer disk metallicities are approximately half those rep- 
resenting the mean value for the inner disk. Therefore, the 
metallicity we measure in extended disks appear to be quan- 
titatively compatible with the predictions of the metallicity 
evolution of the Universe made by recent hydrodynamic sim- 
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ulations. 

It is interesting to compare the metallicity we measure in 
outer disks to the expected metallicit y of Damped Lyman- 
alpha (DLA) systems evolved to z = ( Prochaska et al.|2003[ 
[Kulkarni et al.| '2005). The metallicities of ionized gas in 
emission and neutral gas in absorptio n have been shown to 
agree in at least one local DLA galaxy ( Schulte-Ladbeck et al. 
2005 [l. On face value it appeared that the metal content of 



DLAs fell short of that expected from t he global star forma 
tion rate of galaxies integrated over time fPettini'lQgg*) . How- 
ever, further analysis by Zwaan et al. ( 2005 ) and Chen et al. 
(J2005) found that the metal content ofTbcal field galaxie s Is 
in good agreement with low-z DLAs. Zwaan et al. (2005 i ap- 
plied the luminosity-metallicity relation, together with k nown 
galactic me tallicity gradients (measured only to R25, Vila- 
|Costas &E dmunds 1992), to 21 cm H I maps of local galax- 
ies. They found a median cross section-weighted metallicity 
of [O/H] =-0.85 ± 0.2 (1/7 solar) for H I gas at z = above 



2 X 10^" cm"^). Fig- 



1 in [Zwaan et al. (2005j) shows that typically up to 50% 
he DLA cross-section lies outside R25 . If metallicity gra- 



the DL A column density limit {Nhi 
ure 

of the: 

dients do not continue to decline at the same rate beyond R25 
- as we have found in NGC 1512 and 3621 - a higher me- 
dian would have been measured. Chen et al. (2005) used a 
different approach: they applied average radial H I profiles of 
three galaxy types and measured a Nni-weighted mean metal- 
licity of Z = 0.3 Zq out to the radius where 21 cm emission 
is normally detected. However, not all the gas at these radii 
meets the DLA criteria, so that, had the results been corrected 
for the intrinsic relationship between Nhi and metallicity, a 
higher mean metallicity would have been inferred. 



Utilizing high-quality 2 1 cm maps of NGC 1512 ( Korib- 
alski & Lopez-Sanchez 2009) and NGC 3621 ([Walter et al" 



2008| l, weliave established that every H II region m our sam- 



ple lies in an area satisfying the DLA column density limit. 
We also carried out the same analysis on H II regions in M83 



et al.l 



reported in Bresolin et al. ( 2009b I, overlaying th eir positions 
on the 21 cm maps of Walter et al. (2008i and Crosthwaite 
Ignoring the three Hll regions with R25 > 2.5, 



which are beyond the limit of the H I map, the remaining 46 
H II regions lie in gas above the DLA cutoff. Although our 
sample of three galaxies is small, we are able to compare 
metallicities and H I measure ments in a region of p ara meter 
space not covered directly by |Zwaan et al. ( 2005 1 or Chen 



et al.| d 2005| l. We find that the LO/HJ ratio oTaU gas with 
Nhi = 2x17)20 ^^-2 con-esponds to Z > 0.35 Zq. 

Detailed observations of the metallicity and H I distribution 
beyond R25 need to be performed on a larger sample of galax- 
ies to address this possible mismatch between low redshift 
metal-poor (0.2 solar) DLAs and moderate metallicity (0.35 
solar) outer disks of galaxies. 

7. CONCLUSIONS 

In this paper we have presented the oxygen abundances of a 
large sample of Hll regions in the spiral galaxies NGC 1512 
and NGC 3621, extending out to approximately 25 kpc and 
20 kpc from the galactic centers, respectively. This sample 
allowed us to obtain for the first time the metallicity of the 
extended disks of these two galaxies, and to compare with 
the nebular abundances in the inner disk. Similar work has 
been carried out recently by us and others in a handful of 
other extended-disk galaxies, where the far-UV and Ha emis- 
sion from star-forming regions has been found to large galac- 
tocentric distances. We confirm that a flattening of the ra- 



dial abundances occurs approximately at the isophotal radius, 
which mimicks the behavior of the star formation efficiency. 
The inner disks display the familiar negative metallicity gradi- 
ents, while the outer disks have a virtually constant O/H ratio. 
In our new work, thanks to the large number of high-quality 
spectra obtained, we find what is arguably the 'cleanest' ex- 
ample of such a behavior obtained so far: the abundance scat- 
ter in the extended disk of NGC 3621 is a mere 0.06 dex, 
which is compatible with a perfectly flat metal distribution, 
within the measurement uncertainties. This result does not de- 
pend on the metallicity diagnostic we use. The O/H ratios we 
find in the outer disks of NGC 1512 and NGC 3621 is -0.35 x 
the solar value, based on measurements of the [O III] A4363 
auroral lines, or larger if strong-line metallicity diagnostics 
such as R23 are used. This agrees with our earlier findings 
for M83 and NGC 4625. The statem ent occasionally found in 
the literature (e.g. Espada et al.|201 1 ) that outer disks of spi- 
rals have low metallicities on the order of 0. 1 Zq is therefore 
inaccurate. 

We have summarized a few of the mechanisms that can 
be responsible for the radial abundance distribution and en- 
richment of extended disks. They can be broadly classified 
into the following categories: (a) mixing and turbulence pro- 
cesses; (b) galactic scale outflows and (c) enriched accretion. 
The boundary between internally and externally driven mech- 
anisms is somewhat blurred, because the accretion material is 
progressively enriched by outflows that originate in the galax- 
ies themselves. We can envision a situation in which all these 
processes are at work to produce the chemical signatures we 
observe in extended disks. Moreover, in some cases, e.g. for 
barred or interacting galaxies, the dominant mechanism could 
be different than for structurally different galaxies, such as 
isolated, non-barred galaxies. The two galaxies we have stud- 
ied here are examples of each of these two possibilities, and 
yet the resulting radial abundance distributions (inner slopes, 
outer disk metallicities) are very similar. It is possible that 
a balance between radial mixing and accretion governs the 
chemical evolution of extended disks. 

We have also shown that considering the outer disk metal- 
licities in the context of enriched accretion in cosmologi- 
cal simulations, the empirical values found from H II regions 
can provide useful constraints to the present-time metallicity 
of the infalling gas, in particular the proportion of pristine, 
metal-poor gas relative to the metal enriched gas that origi- 
nates from outflows. More in general, future galactic chem- 
ical evolution models must be able to account for the differ- 
ent radial trends observed in inner and outer disks. In order 
to distinguish between the different scenarios illustrated here, 
specific model predictions for the chemical enrichment and 
the flat gradients in extended disks would be desirable. 

Because of their extreme conditions (e.g. very low gas den- 
sities and long dynamical timescales), outer galaxy disks are 
sensitive probes of the mechanisms leading to the assembly of 
galactic disks and their evolution. The study of the metallic- 
ity of extended disks complements in a unique way the study 
of the star formation properties of outer galaxy disks based 
on UV, optical (Ha) and radio observations, by providing the 
insight into the evolutionary status of galactic systems repre- 
sented by the buildup (and later mixing) of metals in the inter- 
stellar medium. Arguably the best way to measure metallici- 
ties of extended disks to large galactocentric radii beyond the 
nearest galaxies in the foreseeable future is by optical spec- 
troscopy of the very faint H II regions found associated with 
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young stellar clusters. 

There are several questions that are left unanswered regard- 
ing the chemical composition and evolution of extended disks, 
that will be addressed by ongoing and future work. For exam- 
ple, do all galaxies with an extended disk where metallicities 
can be derived from Hll region emission lines behave like 
NGC 1512 and NGC 3621, with an abrupt flattening of the 
abundance gradient roughly at the isophotal radius? Is the ra- 
tio between outer and inner metallicities approximately uni- 
versal? Does it change with galaxy mass? Do the chemical 
properties of outer disks depend on galaxy properties, such 
as mass and environment? Are there parallels in the chem- 



ical mixing properties between extended disks and low sur- 
face brightness galaxies? What is the metal distribution in the 
outer disks of 'rejuvenated' early-type galaxies? 
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APPENDIX 
H II REGION DATA 



In this Appendix we present tables containing the celestial coordinates, the galactocentric distances, line fluxes and extinction 
coefficients for the Hll regions observed in NGC 1512 (Table|5]l and NGC 3621 (Table|6]l. 
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TABLE 5 

NGC 1512: COORDINATES AND REDDENING-CORRECTED FLUXES 



ID 


R.A. 


DEC 


«/«25 


[On] 


[Om] 


[Nil] 


[Sii] 


F(Ha) 




c(H/3) 




(J2000.0) 


(J2000.0) 




3727 


5007 


6583 


6717+6731 


(erg s"' cm 


^) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


1 


04 03 23.9 


-43 16 10.51 


1.94 


260 ± 25 


515 ±25 


20 ±2 


28 ±3 


1.7 xlO" 


5 


0.00 ±0.14 


2 


04 03 52.4 


-43 17 03.36 


1.11 


233 ± 24 


277 ± 14 


21 ±3 


16 ±2 


1.8 xlO" 


5 


0.12 ±0.15 


3 


04 03 25.2 


-43 18 40.37 


1.44 


403 ± 56 


38 ±3 


48 ±8 


57 ±7 


1.7 xlO" 


5 


0.06 ±0.21 


4 


04 04 01.9 


-43 18 59.05 


0.63 


289 ± 28 


47 ±2 


93 ±8 


61 ±4 


1.5 xlO" 


5 


0.32 ±0.14 


5 


04 04 29.6 


-43 19 01.28 


1.58 


363 ± 42 


128 ±7 






1.1 xlO" 


5 


0.00 ±0.17 


6 


04 04 03.1 


-43 19 02.55 


0.64 


275 ± 27 


45 ±2 


91 ±7 


62 ±4 


1.5 xlO" 


5 


0.50 ±0.14 


7 


04 03 35.0 


-43 19 06.02 


1.01 


367 ± 36 


67 ±3 


42 ±5 


54 ±5 


1.0 xlO" 


5 


0.00 ±0.14 


8 


04 04 01.8 


-43 19 10.30 


0.58 


274 ± 30 


41 ±2 


98 ± 12 


70 ±6 


3.1 xlO" 


5 


0.11 ±0.16 


9 


04 04 00.4 


-43 19 11.92 


0.54 


297 ± 25 


88 ±4 


91 ±7 


55 ±3 


6.5 X 10" 


5 


0.44 ±0.12 


10 


04 04 05.7 


-43 19 12.46 


0.67 


185 ± 17 


21 ± 1 


81 ±8 


70 ±5 


4.0 X 10" 


5 


0.07 ±0.13 


11 


04 03 36.2 


-43 19 14.06 


0.94 


300 ± 29 


21 ± 1 


59 ±7 


94 ±8 


1.2 xlO" 


5 


0.01 ± 0.14 


12 


04 03 23.0 


-43 19 16.66 


1.45 


191 ± 15 


515 ±24 


14 ± 1 


18 ± 1 


9.6 X 10- 


5 


0.00 ±0.11 


13 


04 03 56.9 


-43 19 18.21 


0.47 


267 ± 23 


29 ± 1 


105 ±8 


81 ±5 


3.3 X 10" 


5 


0.38 ±0.12 


14 


04 03 57.9 


-43 19 32.45 


0.42 


154 ± 13 


23 ± 1 


97 ±8 


47 ±3 


4.3 X 10" 


5 


0.24 ±0.12 


15 


04 03 38.6 


-43 19 34.33 


0.80 


377 ± 39 


364 ± 18 


44 ±5 


82 ±8 


1.1 xlO" 


5 


0.00 ±0.15 


16 


04 03 53.4 


-43 19 46.93 


0.33 


170 ± 15 


24 ± 1 


123 ±9 


96 ±6 


2.6 X 10" 


5 


0.37 ±0.12 


17 


04 03 52.5 


-43 19 49.85 


0.32 


125 ± 10 


23 ± I 


106 ±7 


46 ±2 


1.2 xlO" 


4 


0.39 ±0.11 


18 


04 04 34.0 


-43 19 56.44 


1.72 


524 ± 57 


168 ±9 






8.0 X 10" 


6 


0.34 ±0.16 


19 


04 03 23.4 


-43 19 58.58 


1.37 


378 ± 39 


231 ± 12 


31 ±3 


38 ±3 


1.3 xlO" 


5 


0.18 ±0.15 


20 


04 03 37.5 


-43 20 05.45 


0.77 


287 ± 25 


73 ±3 


72 ±7 


50 ±4 


2.6 X 10" 


5 


0.04 ±0.12 


21 


04 03 34.4 


-43 20 14.25 


0.89 


421 ± 46 


115±6 


59 ±7 


57 ±5 


6.7 X 10" 


6 


0.11 ±0.16 


22 


04 03 59.9 


-43 20 17.41 


0.29 


87 ±7 


16 ± 1 


106 ± 9 


53 ±3 


1.6 xlO" 


4 


0.24 ±0.12 


23 


04 03 37.2 


-43 20 19.89 


0.76 


255 ± 30 


186 ± 10 


51 ±7 


51 ±5 


1.5 xlO" 


5 


0.01 ± 0.17 


24 


04 04 37.9 


-43 20 19.81 


1.88 


389 ± 36 


176 ± 8 






1.7 XlO" 


5 


0.44 ±0.13 


25 


04 03 50.6 


-43 20 20.79 


0.23 


80 ±7 


7± 1 


1 10 ± 1 1 


52 ±4 


3.5 X 10" 


5 


0.03 ±0.12 


26 


04 03 59.9 


-43 20 26.66 


0.27 


116 ±9 


19 ± 1 


119 ±9 


64 ±4 


9.7 X 10" 


5 


0.30 ±0.11 


27 


04 03 49.0 


-43 20 27.06 


0.27 


134 ± 13 


17 ± 1 


111 ±8 


43 ±2 


1.7 XlO" 


5 


0.56 ±0.13 


28 


04 03 50.0 


-43 20 28.96 


0.23 


155 ± 15 


40 ±2 


108 ±9 


55 ±4 


1.7 XlO" 


5 


0.42 ±0.14 


29 


04 04 37.5 


-43 20 29.22 


1.86 


393 ± 55 


268 ± 15 






2.2 X 10" 


6 


0.36 ±0.21 


30 


04 03 37.8 


-43 20 53.16 


0.70 


313 ±27 


206 ± 10 


89 ±9 


80 ±6 


4.9 X 10" 


5 


0.07 ±0.12 


31 


04 04 34.4 


-43 21 05.13 


1.74 


293 ± 27 


367 ± 17 






1.7 XlO" 


5 


0.36 ±0.13 


32 


04 03 37.6 


-43 21 20.08 


0.71 


319 ± 29 


184 ±9 


61 ±7 


55 ±4 


2.9 X 10" 


5 


0.00 ±0.13 


33 


04 03 37.3 


-43 21 26.05 


0.73 


226 ± 27 


32 ±2 


69 ±9 


65 ±7 


1.1 XlO" 


5 


0.01 ± 0.17 


34 


04 03 33.4 


-43 21 39.03 


0.90 


272 ± 35 


13 ± 1 


75 ± 12 


49 ±7 


3.4 X 10" 


6 


0.00 ±0.19 


35 


04 03 37.7 


-43 21 42.09 


0.73 


361 ± 35 


185 ±9 


61 ±7 


36 ±3 


1.0 XlO" 


5 


0.00 ±0.14 


36 


04 03 47.4 


-43 21 49.09 


0.37 


87 ±7 


17 ± 1 


98 ±9 


57 ±4 


6.0 X 10" 


5 


0.01 ±0.11 


37 


04 03 38.3 


-43 21 49.38 


0.71 


558 ± 59 


74 ±4 


90 ± 12 


97 ± 10 


1.4 XlO" 


5 


0.00 ±0.16 


38 


04 03 50.7 


-43 21 49.60 


0.29 


101 ±9 


14 ± 1 


102 ±8 


49 ±3 


8.8 X 10" 


5 


0.36 ±0.12 


39 


04 03 48.0 


-43 21 56.63 


0.38 


107 ±9 


9± 1 


87 ±7 


53 ±3 


1.3 XlO" 


4 


0.34 ±0.12 


40 


04 03 55.2 


-43 21 59.72 


0.31 


129 ± 10 


24 ± 1 


103 ±9 


67 ±4 


2.9 X 10" 


4 


0.09 ±0.11 


41 


04 03 51.9 


-43 22 03.63 


0.33 


58 ±6 


5± 1 


77 ±8 


47 ±4 


6.1 XlO" 


5 


0.00 ±0.13 


42 


04 03 54.2 


-43 22 05.55 


0.33 


143 ± 12 


22 ± 1 


99 ±9 


69 ±5 


6.5 X 10" 


5 


0.11 ±0.12 


43 


04 03 38.5 


-43 22 06.43 


0.73 


326 ± 28 


78 ±4 


69 ±6 


50 ±3 


4.5 X 10" 


5 


0.22 ±0.12 


44 


04 03 48.2 


-43 22 10.89 


0.43 


159 ± 17 


13 ± 1 


107 ± 1 1 


60 ±5 


7.3 X 10" 


6 


0.28 ±0.15 


45 


04 04 36.3 


-43 22 10.92 


1.87 


248 ± 21 


346 ± 16 






4.8 X 10" 


5 


0.23 ±0.12 


46 


04 03 27.5 


-43 22 18.67 


1.19 




86 ±4 


61 ±6 


73 ±6 


2.3 X 10" 


5 


0.00 ±0.12 


47 


04 03 40.6 


-43 22 21.98 


0.69 


262 ± 26 


51 ±3 


84 ±9 


51 ±4 


3.4 X 10" 


5 


0.11 ±0.14 


48 


04 04 36.9 


-43 22 25.28 


1.91 


368 ± 42 


173 ±9 






1.1 XlO" 


5 


0.32 ±0.17 


49 


04 03 41.3 


-43 22 33.00 


0.70 


435 ± 47 


47 ±3 


84 ±8 


66 ±5 


1.9 XlO" 


5 


0.42 ±0.16 


50 


04 03 41.9 


-43 22 35.83 


0.69 


261 ± 24 


110±5 


82 ±8 


53 ±4 


3.6 X 10" 


5 


0.15 ±0.13 


51 


04 04 37.8 


-43 22 38.24 


1.97 


437 ± 50 


74 ±4 






6.6 X 10" 


6 


0.58 ±0.17 


52 


04 03 52.7 


-43 22 44.27 


0.52 


393 ± 50 


90 ±5 


72 ± 11 


60 ±7 


9.6 X 10" 


6 


0.00 ±0.19 


53 


04 03 46.1 


-43 22 58.66 


0.66 


480 ± 47 


178 ±9 


37 ±4 


65 ±5 


3.0 X 10" 


5 


0.18 ±0.14 


54 


04 03 47.7 


-43 23 11.10 


0.68 


469 ± 39 


113±5 


38 ±4 


63 ±5 


9.9 X 10" 


5 


0.00 ±0.12 


55 


04 03 48.6 


-43 23 13.75 


0.68 


344 ± 20 


231 ± 11 


26 ±2 


35 ±2 


4.6 X 10" 


5 


0.00 ± 0.06 


56 


04 03 50.1 


-43 23 38.99 


0.79 


525 ± 50 


122 ±6 


36 ±4 


35 ±3 


1.7 XlO" 


5 


0.19 ±0.14 


57 


04 03 51.7 


-43 23 42.26 


0.79 


286 ± 23 


291 ± 14 


28 ±3 


38 ±3 


1.5 XlO" 


4 


0.00 ±0.11 


58 


04 03 53.7 


-43 24 08.78 


0.92 


480 ± 67 


160 ±9 


39 ±5 


44 ±4 


2.9 X 10" 


6 


0.40 ± 0.21 


59 


04 03 54.8 


-43 24 16.97 


0.96 


498 ± 57 


90 ±5 


48 ±7 


78 ±8 


5.7 X 10" 


6 


0.00 ±0.17 


60 


04 03 57.6 


-43 24 28.38 


1.03 


553 ± 73 


62 ±4 


41 ±7 


58 ±7 


3.9 X 10" 


6 


0.00 ± 0.20 


61 


04 04 34.0 


-43 24 53.60 


2.12 


266 ± 23 


459 ± 21 






3.0 X 10" 


5 


0.00 ±0.12 


62 


04 04 32.8 


-43 25 21.42 


2.16 


310 ±27 


288 ± 13 






1.5 XlO" 


5 


0.14 ±0.12 



Note. — Units of right ascension are hours, minutes and seconds, and units of dechnation are degrees, arcminutes and arcseconds. 
Line fluxes are in units of H/3= 100. The Ha flux in column (9) has been measured from narrow-band images, and corrected for 
extinction. 
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TABLE 6 

NGC 3621: Coordinates and reddening-corrected fluxes 



ID 


R.A. 


DEC 


«/«25 


[On] 


[Oiii] 


[Nil] 


[Sii] 


F(Ha) 




e(H/3) 




(J2000.0) 


(J2000.0) 




3727 


5007 


6583 


6717+6731 


(erg.s"' cm 


"^) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


1 


1 18 13.8 


-32 39 55.10 


1.97 


261 ± 18 


915 ±43 


18 ± 1 


42 ±3 


4.3 X 10" 


5 


0.00 ± 0.08 


2 


1 18 12.4 


-32 40 06.25 


1.89 


353 ± 19 


192 ±9 


21 ± 1 


45 ±2 


6.6 X 10" 


5 


0.05 ± 0.05 


3 


1 18 08.0 


-32 40 15.60 


1.80 


396 ± 27 


170 ±8 


23 ±2 


36 ±2 


5.5 X 10" 


5 


0.09 ± 0.08 


4 


1 18 05.8 


-32 40 28.99 


1.77 


351 ±28 


280 ± 14 


17 ±2 


36 ±2 


1.5 xlO" 


5 


0.23 ±0.10 


5 


1 18 05.8 


-32 40 44.80 


1.72 


308 ± 34 


495 ± 26 


24 ±3 


63 ±6 


3.9 X 10" 


5 


0.08 ±0.16 


6 


1 18 13.2 


-32 40 51.17 


1.74 


305 ± 36 


86 ±5 


29 ±4 


41 ±5 


6.4 X 10" 


6 


0.00 ±0.17 


7 


1 18 05.5 


-32 40 58.76 


1.68 


206 ± 20 


323 ± 16 


19 ±2 


30 ±3 


2.0 X 10" 


5 


0.19 ±0.13 


8 


1 18 02.5 


-32 42 03.42 


1.59 


321 ± 49 


34 ±3 


42 ±7 


41 ±6 


1.9 xlO" 


6 


0.18 ±0.22 


9 


1 18 02.2 


-32 42 08.81 


1.59 


320 ± 38 


160 ±9 


21 ±4 


31 ± 10 


3.4 X 10" 


(-! 


0.00 ±0.17 


10 


1 17 58.4 


-32 42 19.95 


1.78 


334 ± 18 


184 ±8 


19 ± 1 


29 ± 1 


1.3 xlO" 


4 


0.02 ± 0.05 


11 


1 18 04.9 


-32 42 30.21 


1.43 


353 ± 33 


59 ±3 


35 ±4 


74 ±7 


2.0 X 10" 


5 


0.00 ±0.13 


12 


1 18 02.5 


-32 42 57.15 


1.47 


422 ±31 


176 ±8 


26 ±2 


33 ±2 


2.7 X 10" 


5 


0.27 ± 0.09 


13 


1 18 04.2 


-32 43 34.46 


1.30 


398 ± 36 


91 ±5 


30 ±3 


60 ±5 


2.6 X 10" 


5 


0.00 ±0.12 


14 


1 18 02.1 


-32 44 06.32 


1.37 


551 ±42 


49 ±2 


33 ±2 


63 ±2 


1.2 XlO" 


4 


0.75 ± 0.10 


15 


1 18 08.9 


-32 44 07.42 


1.03 


321 ± 30 


102 ±5 


32 ±4 


56 ±5 


2.6 X 10" 


5 


0.12 ±0.13 


16 


1 17 57.2 


-32 44 17.05 


1.73 


297 ± 15 


231 ± 10 


19 ± 1 


30 ± 1 


6.1 XlO" 


5 


0.13 ±0.04 


17 


1 17 56.3 


-32 44 41.94 


1.80 


345 ± 29 


68 ±3 


25 ±2 


48 ±3 


5.4 X 10" 


5 


0.26 ±0.11 


18 


1 18 03.4 


-32 44 54.01 


1.22 


354 ± 22 


175 ±8 


25 ± 1 


37 ±2 


4.0 X 10" 


5 


0.31 ±0.07 


19 


1 18 13.6 


-32 45 10.91 


0.77 


315 ±23 


57 ±3 


43 ±4 


62 ±4 


3.2 X 10" 


5 


0.00 ± 0.09 


20 


1 18 14.0 


-32 45 13.00 


0.77 


289 ± 14 


251 ± 11 


37 ±2 


69 ±3 


2.3 X 10" 


4 


0.00 ± 0.03 


21 


1 18 01.9 


-32 45 30.09 


1.31 


433 ± 36 


64 ±3 


31 ±3 


54 ±4 


1.1 XlO" 


5 


0.04 ±0.11 


22 


1 18 07.0 


-32 45 49.47 


0.90 


516 ±43 


81±4 


37 ±4 


54 ±4 


2.0 X 10" 


5 


0.00 ±0.11 


23 


1 18 14.1 


-32 46 06.85 


0.57 


247 ± 15 


73 ±3 


68 ±4 


97 ±5 


4.6 X 10" 


4 


0.00 ± 0.06 


24 


1 18 10.5 


-32 46 09.59 


0.65 


469 ± 37 


59 ±3 


46 ±3 


53 ±3 


5.4 X 10" 


5 


0.32 ±0.10 


25 


1 18 07.4 


-32 46 24.25 


0.83 


470 ± 46 


23 ±2 


44±5 


62 ±5 


1.7 XlO" 


5 


0.04 ±0.14 


26 


1 18 07.4 


-32 46 39.53 


0.82 


198 ± 12 


342 ± 16 


27 ± 1 


34 ± 1 


1.8 XlO" 


4 


0.26 ± 0.06 


27 


1 18 07.4 


-32 46 50.79 


0.82 


383 ± 24 


104 ±5 


40 ±3 


54 ±3 


1.1 XlO" 


4 


0.00 ± 0.07 


28 


1 18 09.4 


-32 46 54.84 


0.65 


348 ± 26 


112 ±5 


50 ±4 


53 ±3 


4.8 X 10" 


5 


0.09 ± 0.09 


29 


1 18 14.2 


-32 46 55.88 


0.40 


242 ± 12 


50 ±2 


63 ±3 


59 ±2 


3.1 XlO" 


4 


0.00 ± 0.03 


30 


1 18 14.4 


-32 47 11.69 


0.35 


296 ± 15 


51 ±2 


86 ±4 


60 ±2 


5.1 XlO" 


4 


0.23 ± 0.04 


31 


1 18 09.1 


-32 47 21.04 


0.66 


312 ± 16 


175 ±8 


51 ±3 


84 ±3 


2.8 X 10" 


4 


0.12 ±0.04 


32 


1 18 13.5 


-32 47 31.20 


0.32 


119 ±6 


19 ± 1 


87 ±3 


68 ±2 


2.9 X 10" 


4 


0.41 ± 0.03 


33 


1 18 09.6 


-32 47 39.30 


0.61 


450 ± 48 


134 ±7 


44 ±5 


32 ±3 


5.0 X 10" 


5 


0.00 ±0.15 


34 


1 18 15.6 


-32 47 57.61 


0.19 


185 ± 11 


30 ± 1 


99 ±6 


46 ±2 


2.9 X 10" 


4 


0.00 ± 0.06 


35 


1 18 10.6 


-32 48 02.40 


0.53 


291 ± 14 


113±5 


49 ±2 


43 ± 1 


7.1 XlO" 


4 


0.18 ±0.03 


36 


1 18 12.5 


-32 48 09.63 


0.36 


188 ±9 


40 ±2 


83 ±4 


46 ±2 


5.5 X 10" 


4 


0.00 ± 0.03 


37 


1 18 14.0 


-32 48 33.44 


0.22 


129 ±9 


31 ±2 


95 ±6 


55 ±3 


2.8 X 10" 


4 


0.11 ±0.07 


38 


1 18 18.2 


-32 48 46.25 


0.18 


110±8 


22 ± 1 


100 ±5 


55 ±2 


6.9 X 10" 


4 


0.68 ± 0.08 


39 


1 18 19.0 


-32 49 10.93 


0.22 


149 ±8 


59 ±3 


79 ±3 


68 ±2 


1.4 XlO" 


3 


0.50 ± 0.04 


40 


1 18 24.4 


-32 49 17.86 


0.73 


341 ± 19 


86 ±4 


43 ±3 


51 ±2 


3.5 X 10" 


4 


0.00 ± 0.05 


41 


1 18 21.9 


-32 49 53.13 


0.48 


320 ± 15 


44 ±2 


61 ±3 


45 ± 1 


2.2 X 10" 


3 


0.18 ±0.02 


42 


1 18 23.9 


-32 50 04.77 


0.65 


235 ± 13 


232 ± 11 


39 ±2 


53 ±2 


6.4 X 10" 


4 


0.18 ±0.05 


43 


1 18 12.3 


-32 50 12.93 


0.62 


328 ± 22 


99 ±5 


35 ±3 


33 ±2 


4.3 X 10" 


4 


0.00 ± 0.08 


44 


1 18 25.5 


-32 50 26.18 


0.80 


386 ± 23 


83 ±4 


44 ±3 


41 ±2 


3.9 X 10" 


4 


0.00 ± 0.06 


45 


1 18 24.0 


-32 5102.15 


0.69 


425 ± 45 


62 ±4 


45 ±5 


48 ±4 


1.2 XlO" 


4 


0.00 ±0.15 


46 


1 18 22.2 


-32 5139.81 


0.65 


162 ±8 


399 ± 18 


18 ± 1 


17 ± 1 


5.0 X 10" 


4 


0.00 ± 0.02 


47 


1 18 19.7 


-32 5145.83 


0.61 


413 ±39 


28 ±2 


46 ±5 


65 ±5 


1.1 XlO" 


4 


0.00 ±0.13 


48 


1 18 20.4 


-32 52 28.62 


0.76 


366 ± 23 


61 ±3 


43 ±3 


44 ±2 


6.3 X 10" 


5 


0.20 ± 0.07 


49 


1 18 32.8 


-32 52 31.07 


1.46 


276 ± 14 


340 ± 15 


15 ± 1 


25 ± 1 


1.3 XlO" 


4 


0.13 ±0.04 


50 


1 18 32.2 


-32 52 44.55 


1.42 


350 ± 24 


278 ± 13 


24 ±2 


44 ±2 


9.4 X 10" 


5 


0.07 ± 0.08 


51 


1 18 28.9 


-32 52 51.41 


1.17 


166 ±8 


357 ± 16 


15 ± 1 


22 ± 1 


2.0 X 10" 


4 


0.15 ±0.03 


52 


1 18 32.8 


-32 53 01.62 


1.48 


241 ± 12 


310 ± 14 


14 ± 1 


22 ± 1 


4.2 X 10" 


4 


0.02 ± 0.03 


53 


1 18 28.8 


-32 53 07.35 


1.19 


160 ±7 


445 ± 20 


11 ± 1 


17 ±1 


2.4 X 10" 


3 


0.14 ±0.01 


54 


1 18 33.2 


-32 53 49.27 


1.55 


255 ± 19 


293 ± 14 


18 ± 1 


47 ±3 


4.9 X 10" 


5 


0.24 ± 0.09 


55 


1 18 28.9 


-32 53 49.29 


1.27 


419 ±38 


33 ±2 


27 ±2 


62 ±4 


2.3 X 10" 


5 


0.25 ±0.12 


56 


1 18 28.6 


-32 54 04.33 


1.29 


284 ± 17 


152 ±7 


23 ± 1 


37 ±2 


1.3 XlO" 


4 


0.09 ± 0.06 


57 


1 18 28.7 


-32 54 15.30 


1.31 


406 ± 29 


103 ±5 


29 ±2 


95 ±5 


7.5 X 10" 


5 


0.17 ±0.08 


58 


1 18 33.1 


-32 54 21.71 


1.59 


370 ± 38 


164 ±8 


16 ±2 


37 ±3 


8.8 XlO" 


5 


0.21 ±0.15 


59 


1 18 28.3 


-32 54 22.05 


1.31 


365 ± 27 


166 ±8 


25 ±1 


45 ±2 


5.4 X 10" 


5 


0.68 ± 0.09 


60 


1 18 27.8 


-32 54 30.96 


1.31 


282 ± 16 


398 ± 18 


21 ±1 


46 ±1 


1.5 XlO" 


4 


0.53 ± 0.05 



Note. — (Continues) 
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TABLE 6 
NGC 3621: Coordinates and reddening-corrected fluxes 



ID 


R.A. 


DEC 


«/«25 


[On] 


[Oiii] 


[Nil] 


[Sii] 


F(Hq) 




c(H/3) 




(J2000.0) 


(J2000.0) 




3727 


5007 


6583 


6717+6731 


(erg s"' cm 


^) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


61 


1 18 18.5 


-32 54 31.71 


1.25 


428 ± 22 


102 ±5 


31 ± 1 


67 ±2 


2.1 xlO" 


14 


0.38 ± 0.04 


62 


1 18 27.8 


-32 54 37.12 


1.33 


317 ± 17 


268 ± 12 


23 ± 1 


40 ± 1 


2.5 X 10" 


14 


0.36 ± 0.04 


63 


1 18 27.8 


-32 54 46.24 


1.35 


294 ± 15 


372 ± 17 


16 ± 1 


30 ± 1 


1.9 xlO" 




0.48 ± 0.04 


64 


1 18 23.4 


-32 55 06.61 


1.32 


273 ± 13 


352 ± 16 


16 ± 1 


27 ± 1 


3.0 X 10" 




0.05 ± 0.03 


65 


1 18 26.6 


-32 55 11.70 


1.38 


403 ± 52 


89 ±4 


23 ±2 


40 ±3 


8.4 X 10" 




0.60 ± 0.20 


66 


1 18 23.5 


-32 55 21.06 


1.37 


197 ± 15 


446 ± 21 


14 ± 1 


32 ±2 


4.3 X 10" 




0.17 ±0.09 


67 


1 18 24.8 


-32 55 35.23 


1.43 


383 ± 21 


214 ± 10 


22 ± 1 


47 ±2 


6.5 X 10" 




0.37 ± 0.05 


68 


1 18 15.1 


-32 55 50.76 


1.69 


411 ±37 


118±6 


21 ±2 


34 ±3 


3.8 X 10" 




0.14 ±0.12 


69 


1 18 24.6 


-32 55 54.09 


1.49 


287 ± 30 


131 ±7 


27 ±3 


44 ±3 


1.0 xlO" 




0.33 ±0.15 


70 


1 18 23.1 


-32 56 10.95 


1.54 


441 ± 39 


37 ±2 


27 ±2 


51 ±3 


3.9 X 10" 




0.34 ±0.12 


71 


1 18 20.6 


-32 56 34.97 


1.66 


581 ±60 


103 ±5 


19 ±2 


43 ±3 


2.2 X 10" 




0.45 ±0.15 


72 


1 18 20.0 


-32 57 08.01 


1.81 


455 ± 27 


75 ±3 


22 ± 1 


50 ±2 


8.2 X 10" 




0.11 ±0.06 


73 


1 18 17.5 


-32 57 24.66 


1.95 


452 ± 33 


133 ±6 


22 ±2 


60 ±4 


4.5 X 10" 




0.00 ± 0.09 



Note. — Units of right ascension are hours, minutes and seconds, and units of dechnation are degrees, arcminutes and arcseconds. 
Line fluxes are in units of H/3 = 100. The Ha flux in column (9) has been measured from narrow -band images, and corrected for 

extinction. 



